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ABSTRACT 

-lv 

^>"This  paper  presents  approximation  theory  for  the  linear-quadratic- 
Gaussian  optimal  control  problem  for  flexible  structures  whose  distributed 
models  have  bounded  input  and  output  operators.  The  main  purpose  of  the 
theory  is  to  guide  the  design  of  finite  dimensional  compensators  that  approxi¬ 
mate  closely  the  optimal  compensator  separates  into  an  optimal  linear- 
quadratic  control  problem  lies  in  the  solution  to  an  infinite  dimensional 
Riccati  operator  equation.  The  approximation  scheme  in  the  paper  approximates 
the  infinite  dimensional  LQG  problem  with  a  sequence  of  finite  dimensional  LQG 
problems  defined  for  a  sequence  of  finite  dimensional,  usually  finite  element 
or  modal,  approximations  of  the  distributed  model  of  the  structure.  Two 
Riccati  matrix  equations  determine  the  solution  to  each  approximating  problem. 

The  finite  dimensional  equations  for  numerical  approximation  are  de¬ 
veloped,  including  formulas  for  converting  matrix  control  and  estimator  gains 
to  their  functional  representation  to  allow  comparison  of  gains  based  on  dif¬ 
ferent  orders  of  approximation.  Convergence  of  the  approximating  control  and 
estimator  gains  and  of  the  corresponding  finite  dimensional  compensators  is 
studied.  Also,  convergence  and  stability  of  the  closed-loop  systems  produced 
with  the  finite  dimensional  compensators  are  discussed,  v  The  convergence 
theory  is  based  on  the  convergence  of  the  solutions  of  the  finite  dimensional 
Riccati  equations  to  the  solutions  of  the  infinite  dimensional  Riccati  equa¬ 
tions.  A  numerical  example  with  a  flexible  beam,  a  rotating  rigid  body,  and  a 
lumped  mass  is  given. 
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1.  Introduction 


Recent  years  have  seen  increasing  research  in  active  control  of  flexible 
structures.  The  primary  motivation  for  this  research  is  oontrol  of  large 
flexible  aerospace  structures,  which  are  becoming  larger  and  more  flexible  at 
the  same  time  that  their  performance  requirements  are  becoming  more  stringent. 
For  example,  in  tracking  and  other  applications,  satellites  with  large  anten¬ 
nae,  solar  collectors  and  other  flexible  components  must  perform  fast  slew 
maneuvers  while  maintaining  tight  control  over  the  vibrations  of  their  flexi¬ 
ble  elements.  Both  of  these  conflicting  objectives  can  be  achieved  only  with 
a  sophisticated  controller.  There  are  applications  also  to  control  of  robotic 
manipulators  with  flexible  links,  and  possibly  to  stabilization  of  large  civil 
engineering  structures  such  as  long  bridges  and  tall  buildings. 

The  first  question  that  must  be  answered  when  designing  a  controller  for 
a  flexible  structure  is  whether  a  finite  dimensional  model  is  sufficient  as  a 
basis  for  a  controller  that  will  produce  the  required  performance,  or  is  a 
distributed  model  necessary?  While  some  structures  can  be  modeled  well  by  a 
fixed  number  of  dominant  modes,  there  are  structures  whose  flexible  character 
can  be  captured  sufficiently  for  precise  control  only  by  a  distributed  model. 
Still  others  —  perhaps  most  of  the  aerospace  structures  of  the  future  —  can 
be  modeled  sufficiently  for  control  purposes  by  some  finite  dimensional 
approximation,  but  an  adequate  approximation  may  be  impossible  to  determine 
before  design  of  the  controller,  or  compensator.  This  paper  deals  with  struc¬ 
tures  that  are  flexible  enough  to  require  a  distributed  model  in  the  design  of 


an  optimal  LQG  compensator. 


The  linear-quadratic-gaussian  optimal  control  problem  for  distributed,  or 
infinite  dimensional,  systems  is  a  generalization  to  Hilbert  space  of  the  now 
classical  LQG  problem  for  finite  dimensional  systems.  The  solution  to  the 
infinite  dimensional  problem  yields  an  infinite  dimensional  state- estimator- 
based  compensator,  which  is  optimal  in  the  context  of  this  paper.  By  a 
separation  principle  [Bl,  C4]  ,  the  problem  reduces  to  a  deterministic  linear- 
quadratic  optimal  control  problem  and  an  optimal  estimation,  or  filtering, 
problem  with  gaussian  white  noise.  In  infinite  dimensions,  the  control  system 
dynamics  are  represented  by  a  semigroup  of  bounded  linear  operators  instead  of 
the  matrix  exponential  operators  in  finite  dimensions,  and  the  plant  noise 
process  may  be  an  infinite  dimensional  random  process.  The  solutions  to  both 
the  control  and  filtering  problems  involve  Riccati  operator  equations,  which 
are  generalizations  of  the  Riccati  matrix  equations  in  the  finite  dimensional 
case.  Current  results  on  the  infinite  dimensional  LQG  problem  are  most  com¬ 
plete  for  problems  where  the  input  and  measurement  operators  are  bounded,  as 
this  paper  requires  throughout.  This  boundedness  also  permits  the  strongest 
approximation  results  here.  For  related  control  problems  with  unbounded  input 
and  measurement,  see  [C3,  C5 ,  LI,  L2]  . 

Our  primary  objective  in  this  paper  is  to  approximate  the  optimal  infin¬ 
ite  dimensional  LQG  compensator  for  a  distributed  model  of  a  flexible  struc¬ 
ture  with  finite  dimensional  compensators  based  on  approximations  to  the 
structure,  and  to  have  these  finite  dimensional  compensators  produce  near 
optimal  performance  of  the  closed-loop  system.  We  discuss  how  the  gains  that 
determine  the  finite  dimensional  compensators  converge  to  the  gains  that 
determine  the  infinite  dimensional  compensator,  and  we  examine  the  sense  in 
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which  the  finite  dimensional  compensators  converge  to  the  infinite  dimensional 
compensator.  With  this  analysis,  we  can  predict  the  performance  of  the 
closed-loop  system  consisting  of  the  distributed  plant  and  a  finite  dimen¬ 
sional  compensator  that  approximates  the  infinite  dimensional  compensator. 

Our  design  philosophy  is  to  let  the  convergence  of  the  finite  dimensional 
compensators  indicate  the  order  of  the  compensator  that  is  required  to  produce 
the  desired  performance  of  the  structure.  The  two  main  factors  that  govern 
rate  of  convergence  are  the  desired  performance  (e.g. ,  fast  response)  and  the 
structural  damping.  We  should  note  that  any  one  of  our  compensators  whose 
order  is  not  sufficient  to  approximate  the  infinite  dimensional  compensator 
closely  will  not  in  general  be  the  optimal  compensator  of  that  fixed  order; 
i.e. ,  the  optimal  fixed-order  compensator  that  would  be  constructed  with  the 
design  philosophy  in  [B7 ,  B8] .  But  as  we  increase  the  order  of  approximation 
to  obtain  convergence,  our  finite  dimensional  compensators  become  essentially 
identical  to  the  compensator  that  is  optimal  over  compensators  of  all  orders. 

An  important  question,  of  course,  is  how  large  a  finite  dimensional  com¬ 
pensator  we  must  use  to  approximate  the  infinite  dimensional  compensator.  In 
[G6,  G7,  G8 ,  Ml],  we  have  found  that  our  complete  design  strategy  yields  com¬ 
pensators  of  reasonable  size  for  distributed  models  of  complex  space  struc¬ 
tures.  This  strategy  in  general  requires  two  steps  to  obtain  an  implementable 
compensator  that  is  essentially  identical  to  the  optimal  infinite  dimensional 
compensator:  the  first  step  determines  the  optimal  compensator  by  letting  the 
finite  dimensional  compensators  converge  to  it;  the  second  step  reduces,  if 
possible,  the  order  of  a  large  (converged)  approximation  to  the  optimal  com¬ 
pensator.  The  first  step,  which  is  the  one  involving  control  theory  and 
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approximation  theory  for  distributed  systems,  is  the  subject  of  this  paper. 

For  the  second  step,  a  simple  modal  truncation  of  the  large  compensator  some¬ 
times  is  sufficient,  but  there  are  more  sophisticated  methods  in  finite  dimen¬ 
sional  control  theory  for  order  reduction.  For  example  [G8,  Ml],  we  have 
found  balanced  realizations  [M2]  to  work  well  for  reducing  large  compensators. 

The  approximation  theory  in  this  paper  follows  from  the  application  of 
approximation  results  in  [B6,  G3 ,  G4]  to  a  sequence  of  finite  dimensional 
optimal  LQG  problems  based  on  a  Ritz-Galerkin  approximation  of  the  flexible 
structure.  For  the  optimal  linear-quadratic  control  problem,  the  approxima¬ 
tion  theory  here  is  a  substantial  improvement  over  that  in  [Gl]  because  here 
we  allow  rigid-body  modes,  more  general  structural  damping  (including  damping 
in  the  boundary),  and  much  more  general  finite  element  approximations.  These 
generalizations  are  necessary  to  accommodate  common  features  of  oomplex  space 
structures  and  the  most  useful  finite  element  schemes.  For  example,  we  write 
the  equations  for  constructing  the  approximating  control  and  estimator  gains 
and  finite  dimensional  compensators  in  terms  of  matrices  that  are  built 
directly  from  typical  mass,  stiffness  and  damping  matrices  for  flexible  struc¬ 
tures,  along  with  actuator  influence  matrices  and  measurement  matrices. 

For  the  estimator  problem,  this  paper  presents  the  first  rigorous  approx¬ 
imation  theory.  (We  have  used  less  complete  versions  of  the  results  in  previ¬ 
ous  research  [G6 ,  G7 ,  G8 ,  Ml].)  As  in  the  finite  dimensional  case,  the  infin¬ 
ite  dimensional  optimal  estimation  problem  is  the  dual  of  the  infinite  dimen¬ 
sional  optimal  control  problem,  and  the  solutions  to  both  problems  have  the 
same  structure.  Because  we  exploit  this  duality  to  obtain  the  approximation 
theory  for  the  estimation  problem  from  the  approximation  theory  for  the 
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optimal  control  problem,  the  analysis  in  this  paper  is  almost  entirely  deter¬ 
ministic.  We  discuss  the  stochastic  interpretation  of  the  estimation  problem 
and  the  approximating  state  estimators  briefly,  but  we  are  concerned  mainl> 
with  deterministic  questions  about  the  structure  and  convergence  of  approxima¬ 
tions  to  an  infinite  dimensional  compensator  and  the  performance  —  especially 
stability  —  of  the  closed-loop  systems  produced  by  the  approximating  compen¬ 
sators. 

Next,  an  outline  of  the  paper  should  help.  The  paper  has  two  main  parts, 
which  correspond  roughly  to  the  separation  of  the  optimal  LAG  problem  into  an 
optimal  linear-quadratic  regulator  problem  and  an  optimal  state  estimation 
problem.  The  first  half.  Sections  2  through  6,  deal  with  the  control  system 
and  the  optimal  regulator  problem.  Sections  7  through  10  treat  the  state 
estimator  and  the  compensator  that  is  formed  by  applying  the  control  law  of 
the  first  half  of  the  paper  to  the  output  of  the  estimator. 

Section  2  defines  the  abstract  model  of  a  flexible  structure  and  the 
energy  spaces  to  be  used  throughout  the  paper.  We  assume  a  finite  number  of 
actuators,  since  this  is  the  case  in  all  applications,  and  we  assume  that  the 
actuator  influence  operator  is  bounded.  Section  2  also  establishes  certain 
mathematical  properties  of  the  open-loop  system  that  are  useful  in  control  and 
approximation.  To  our  knowledge,  the  exponential  stability  theorem  in  Section 
2.3  is  a  new  result,  and  we  find  it  interesting  that  such  a  simple  Lyapunov 
functional  accommodates  such  a  general  damping  model. 

Section  3  discusses  the  linear-quadratic  optimal  control  problem  for  the 
distributed  model  of  the  structure  and  establishes  some  estimates  involving 
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bounds  on  solutions  to  infinite  dimensional  Riccati  equations  and  open-loop 
and  closed-loop  decay  rates.  We  need  these  estimates  for  the  subsequent 
approximation  theory.  To  get  the  approximation  theory  for  the  estimation 
problem,  we  have  to  give  certain  results  on  the  control  problem  in  a  more  gen¬ 
eral  form  than  would  be  necessary  were  we  interested  only  in  the  control  prob¬ 
lem  for  flexible  structures.  Therefore,  in  Section  3,  as  in  Sections  5  and  7, 
we  first  give  some  generic  results  applicable  to  the  LQG  problem  for  a  variety 
of  distributed  systems  and  then  apply  the  generic  results  to  the  control  of 
flexible  structures. 

Because  we  assume  a  finite  number  of  actuators  and  a  bounded  input  opera¬ 
tor,  the  optimal  feedback  control  law  consists  of  a  finite  number  of  bounded 
linear  functionals  on  the  state  space,  which  is  a  Hilbert  space.  This  means 
that  that  the  feedback  law  can  be  represented  in  terms  of  a  finite  number  of 
vectors,  which  we  call  functional  control  gains,  whose  inner  products  with  the 
generalized  displacement  and  velocity  vectors  define  the  control  law.  For  any 
finite-rank,  bounded  linear  feedback  law  for  a  control  system  on  a  Hilbert 
space,  the  existence  of  3uch  gains  is  obvious  and  well  known.  A  functional 
control  gain  for  a  flexible  structure  will  have  one  or  more  distributed  com¬ 
ponents,  or  kernels,  corresponding  to  each  distributed  component  of  the  struc¬ 
ture  and  scalar  components  corresponding  to  each  rigid  component  of  the  struc¬ 
ture. 


We  introduce  the  functional  control  gains  at  the  end  of  Section  3,  and  we 
introduce  analogous  functional  estimator  gains  in  Section  7.  The  functional 
gains  play  a  prominent  roll  in  our  analysis.  They  give  a  concrete  representa¬ 
tion  of  the  infinite  dimensional  compensator  and  provide  a  criterion  for  con- 
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vergenee  of  the  approximating  finite  dimensional  compensators. 

We  develop  the  approximation  scheme  for  the  control  problem  in  Section  4. 
The  idea  is  to  solve  a  finite  dimensional  linear-quadratic  regulator  problem 
for  each  of  a  sequence  of  Ritz-Galerkin  approximations  to  the  structure.  We 
develop  the  approximation  of  the  structure  in  Section  4.1  and  prove  conver¬ 
gence  of  the  approximating  open-loop  systems.  The  approximation  scheme 
includes  most  finite  element  approximations  of  flexible  structures.  For  con¬ 
vergence,  we  use  the  Trotter-Kato  semigroup  approximation  theorem,  which  was 
used  in  optimal  open-loop  control  of  hereditary  systems  in  [B5]  and  has  been 
used  in  optimal  feedback  control  of  hereditary,  hyperbolic  and  parabolic  sys¬ 
tems  in  [B6,  G1 ,  G3]  and  other  papers.  The  usual  way  to  invoke  Trotter-Kato 
is  to  prove  that  the  resolvents  of  the  approximating  semigroup  generators  con¬ 
verge  strongly.  To  prove  this,  we  introduce  an  inner  product  that  involves 
both  the  strain-energy  inner  product  and  the  damping  functional,  and  show  that 
the  resolvent  of  each  finite  dimensional  semigroup  generator  is  the  projec¬ 
tion,  with  respect  to  this  special  inner  product,  of  the  resolvent  of  the  ori¬ 
ginal  semigroup  generator  onto  the  approximation  subspace.  The  idea  works  as 
well  for  the  adjoints  of  the  resolvents,  and  when  the  open-loop  semigroup  gen¬ 
erator  has  compact  resolvent,  it  follows  from  our  projection  that  the  approxi¬ 
mating  resolvent  operators  converge  in  norm. 

In  Section  4.2,  we  define  the  sequence  of  finite  dimensional  optimal  con¬ 
trol  problems,  whose  solutions  approximate  the  solution  to  the  infinite  dimen¬ 
sional  problem  of  Section  3.  The  solution  to  each  finite  dimensional  problem 
is  based  on  the  solution  to  a  Riccati  matrix  equation,  and  we  give  formulas 
for  using  the  solution  to  the  Riccati  matrix  equation  to  compute  approxima- 
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tions  to  the  functional  control  gains  as  linear  combinations  of  the  basis  vec¬ 
tors. 


Section  5.1  summarizes  some  generic  convergence  results  from  [B6,  G3  ,  G4] 
on  approximation  of  solutions  to  infinite  dimensional  Riccati  equations.  Sec¬ 
tion  5.2  applies  these  generic  results  to  obtain  sufficient  conditions  for 
convergence  —  and  nonconvergence  —  of  the  solutions  of  the  approximating 
optimal  control  problems  in  Section  4.2.  A  main  sufficient  condition  for  con¬ 
vergence  is  that  the  structure  have  damping,  however  small,  that  makes  all 
elastic  vibrations  of  the  open-loop  system  exponentially  stable.  This  is  a 
necessary  condition  if  the  state  weighting  operator  in  the  control  problem  is 
coercive. 

In  Section  6,  we  present  an  example  in  which  the  structure  consists  of  an 
Euler^- Bernoulli  beam  attached  on  one  end  to  a  rotating  rigid  hub  and  on  the 
other  end  to  a  lumped  mass.  We  emphasize  the  fact  that  we  do  not  solve,  or 
even  write  down,  the  coupled  partial  and  ordinary  differential  equations  of 
motion.  For  both  the  definition  and  numerical  solution  of  the  problem,  only 
the  kinetic  and  strain  energy  functionals  and  a  dissipation  functional  for  the 
damping  are  required.  We  show  the  approximating  functional  control  gains 
obtained  by  using  a  standard  finite  element  approximation  of  the  beam,  and  we 
discuss  the  effect  on  convergence  of  structural  damping  and  of  the  ratio  of 
state  weighting  to  control  weighting  in  the  performance  index.  As  suggested 
by  a  theorem  in  Section  5,  the  functional  gains  do  not  converge  when  no  struc¬ 
tural  damping  is  modeled. 

In  Section  7,  we  begin  the  theory  for  closing  the  loop  on  the  control 
system.  We  assume  a  finite  number  of  bounded  linear  measurements  and 
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construct  the  optimal  state  estimator,  which  is  infinite  dimensional  in  gen¬ 
eral.  The  gains  for  this  estimator  are  obtained  from  the  solution  to  an 
infinite  dimensional  Riccati  equation  that  has  the  same  form  as  the  infinite 
dimensional  Riccati  equation  in  the  control  problem.  We  call  these  gains 
functional  estimator  gains  because  they  are  vectors  in  the  state  space. 

Since  the  approximation  issues  that  this  paper  treats  are  fundamentally 
deterministic,  we  make  the  paper  self  contained  by  defining  the  infinite 
dimensional  estimator  as  an  observer,  although  the  only  justification  for  cal¬ 
ling  this  estimator  and  the  corresponding  compensator  optimal  is  their 
interpretation  in  the  context  of  stochastic  estimation  and  control.  We  dis¬ 
cuss  the  stochastic  interpretation  but  do  not  use  it.  We  say  estimator  and 
observer  interchangeably  to  emphasize  the  deterministic  definition  of  the 
estimator  here. 

With  the  optimal  control  law  of  Section  3  and  the  optimal  estimator  of 
Section  7,  we  construct  the  optimal  compensator,  which  also  is  infinite  dimen¬ 
sional  in  general.  The  transfer  function  of  this  compensator  is  irrational, 
but  it  is  still  an  m(number  of  actuators)  p (number  of  sensors)  matrix  func¬ 
tion  of  a  complex  variable,  as  in  finite  dimensional  control  theory.  The 
optimal  closed-loop  system  consists  of  the  distributed  model  of  the  structure 
controlled  by  the  optimal  compensator. 

Approximation  of  the  optimal  compensator  is  based  on  approximating  the 
infinite  dimensional  estimator  with  the  sequence  of  finite  dimensional  estima¬ 
tors  defined  in  Section  8.1.  The  gains  for  the  approximating  estimators  are 
given  in  terms  of  the  solutions  to  finite  dimensional  Riccati  equations  that 
approximate  the  infinite  dimensional  Riccati  equation  in  Section  7.  Although 
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defined  as  observers,  these  finite  dimensional  estimators  can  be  interpreted 
as  Kalman  filters,  as  shown  in  Section  8.2.  In  Section  8.3,  we  give  formulas 
for  finite  dimensional  functional  estimator  gains  that  approximate  the  func¬ 
tional  estimator  gains  of  Section  7.  These  approximating  estimator  gains 
indicate  how  closely  the  finite  dimensional  estimators  approximate  the  infin¬ 
ite  dimensional  estimator.  In  Section  8.4,  we  apply  the  Riccati  equation 
approximation  theory  of  Section  5  to  describe  the  convergence  of  the  finite 
dimensional  estimators. 

Most  of  the  results  in  Section  8  are  analogous  to  results  for  the  control 
problem  and  follow  from  the  same  basic  approximation  theory,  but  certain 
differences  require  careful  analysis.  There  is  an  important  difference  in  the 
way  that  the  Riccati  matrices  to  be  computed  are  defined  in  terms  of  the  fin¬ 
ite  dimensional  Riccati  operators.  Indeed,  the  Riccati  matrix  equations  to  be 
solved  numerically  might  seen  incorrect  at  first.  To  demonstrate  that  the 
finite  dimensional  estimators  that  we  define  in  Section  8.1  are  natural 
approximations  to  the  optimal  infinite  dimensional  estimator,  we  show  in  Sec¬ 
tion  8.2  that  each  finite  dimensional  estimator  is  a  Kalman  filter  for  the 
corresponding  finite  element  approximation  of  the  flexible  structure.  The 
brief  discussion  in  Section  8.2  is  the  only  place  in  the  paper  where  stochas¬ 
tic  estimation  theory  is  necessary,  and  none  of  the  analysis  in  the  rest  of 
the  paper  depends  on  this  discussion. 

In  Section  9.1,  we  apply  the  nth  control  law  of  Sect: on  4  to  the  output 
of  the  nth  estimator  to  form  the  nth  compensator.  (The  order  of  approximation 
is  n. )  The  nth  closed-loop  system  consists  of  the  distributed  model  of  the 
structure  controlled  by  the  nth  compensator.  Since  each  finite  dimensional 
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estimator  is  realizable,  the  nth  compensator  and  the  nth  closed-loop  system 


are  realizable.  In  Section  9.2,  we  discuss  how  the  sequence  of  realizable 
closed-loop  systems  approximates  the  optimal  closed-loop  system.  Probably  the 
most  important  question  here  is  whether  exponential  stability  of  the  optimal 
closed-loop  system  implies  exponential  stability  of  the  nth  closed-loop  system 
for  n  sufficiently  large.  We  have  been  able  to  prove  this  only  when  the 
approximation  basis  vectors  are  the  natural  modes  of  undamped  free  vibration 
and  these  modes  are  not  coupled  by  structural  damping.  That  this  stability 
result  can  be  generalized  is  suggested  by  the  results  in  Section  9.3,  which 
describe  how  the  transfer  functions  of  the  finite  dimensional  compensators 
approximate  the  transfer  function  of  the  optimal  compensator. 

In  Section  10,  we  complete  the  compensator  design  for  the  example  in  Sec¬ 
tion  6.  Assuming  that  white  noise  corrupts  the  single  measurement  and  that 
distributed  white  noise  disturbs  the  structure,  we  compute  the  gains  for  the 
finite  dimensional  estimators  and  show  the  functional  estimator  gains.  As  in 
the  control  problem,  the  functional  gains  do  not  converge  when  no  damping  is 
modeled.  We  apply  the  control  laws  computed  in  Section  6  to  the  output  of  the 
estimators  in  Section  10  to  construct  the  finite  dimensional  compensators,  and 
we  show  the  frequency  response  of  these  compensators.  As  predicted  by  Section 
9.3,  the  frequency  response  of  the  nth  compensator  converges  to  the  frequency 
response  of  the  optimal  infinite  dimensional  compensator  as  n  increases.  In 
Section  10.3,  we  discuss  the  structure  and  dimension  of  the  finite  dimensional 
compensator  that  should  be  implemented. 

We  conclude  in  Section  11  by  discussing  where  the  approximation  theory 
presented  in  this  paper  is  most  complete  and  what  further  results  would  be 
most  important. 
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2.  The  Control  System 

We  consider  the  system 

x(t>  +  D0x(t)  +  AQx(t) 
where  x(t)  is  in  a  real  Hilbert  space 

The  linear  stiffness  operator  is  densely  defined  and  selfadjoint  with  com¬ 
pact  resolvent  and  at  most  a  finite  number  of  negative  eigenvalues.  We  will 
postpone  discussion  of  the  damping  operator  DQ  momentarily,  except  to  say  that 
it  is  symmetric  and  nonnegative.  The  input  operator  BQ  is  a  linear  operator 
from  R®  to  H,  hence  bounded. 

By  natural  modes,  m§  will  mean  J&s  eigenvectors  ^  <?f  eigenvalue 


*  BqU ( t ) ,  t  >  0, 


(2.1) 


H  and  u(t)  is  in  R™  for  some  finite  m. 


Vj  '  Vj  * 


(2.2) 


From  our  hypotheses  on  Aq,  we  know  that  these  eigenvalues  form  an  infinitely 
increasing  sequence  of  real  numbers,  of  which  all  but  a  finite  number  are 
positive.  Also,  the  corresponding  eigenvectors  are  complete  in  H  and  satisfy 


<^i*  “  <A0^i'  1  * 


(2.3) 


(These  properties  of  the  eigenvalue  problem  (2.2)  are  standard.  See,  for 
example,  [Bl]  ,  [Kll  .)  For  Aj  >  0,  u>j  =  “f  ^  is  a  natural  frequency. 
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Remark  2*1*  Our  analysis  Includes  the  system 


H»x(t)  +  Dnx(t)  +  JLx(t)  =  Bnu(t),  t  >  0, 

u  u  u  u  (2.1  ; 

where  the  mass  operator  Mq  is  a  selfadjoint,  bounded  and  coercive  linear 

operator  on  a  real  Hilbert  space  HQ.  The  operators  A0,  Bq  and  DQ  in  (2.1') 

have  the  same  properties  with  respect  to  Hq  that  the  corresponding  operators 

in  (2.1)  have  with  respect  to  H.  To  include  (2.1')  in  our  analysis,  ve  need 

only  take  H  to  be  HQ  with  the  norm- equivalent  inner  product  = 

<M0\‘>h  ,  and  multiply  (2.1')  on  the  left  by  Mq*.  In  H,  the  operator  M^Aq 

is  selfadjoint  with  compact  resolvent,  and  Mq*Dq  is  symmetric  and  nonnegative. 
With  no  loss  of  generality  then,  we  will  refer  henceforth  only  to  (2.1)  and 
assume  that  the  H-inner  product  accounts  for  the  the  mass  distribution.  Q 


2.1  The  Energy  Spaces  and  the  First-order  Form  of  the  System 


IllS  Elastic-Strain- Energy  Space  V  and  Total-Energy  Space  E 


A0  + 


We  choose  a  bounded,  selfadjoint  linear  operator  Aj  on  H  such  that  Aq  = 
A^  is  ooercive;  i.e. ,  there  exists  p  >  0  for  which 

<Aqx,x>h  2  p||x||*,  x  e  D(Aq)  =  D(Aq) .  ^ 


Since  Aq  is  bounded  from  below,  there  will  be  infinitely  many  such  A^'s. 
In  applications  like  our  example  in  Section  6,  it  is  natural  to  select  for  Aj 
an  operator  whose  null  space  is  the  orthogonal  complement  (in  H)  of  the  eigen- 
space  of  Aq  corresponding  to  nonpositive  eigenvalues.  Obviously,  any  that 
makes  AQ  coercive  must  be  positive  definite  on  the  nonpositive  eigenspaoe  of 


With  chosen,  we  define  the  Hilbert  space  V  to  be  the  completion  of 
D(Ag)  with  respect  to  the  inner  product  =  <AqVj,V2>jj»  v^  and 

v2  c  D(A0).  Note  that  V  =  D(Aj^2)  and  <vi>v2^y  =  *^0^2v2>H*  ^Sinee  Ai 

is  a  bounded  operator  on  H,  different  choices  of  Aj  yield  V's  with  equivalent 
norms,  thus  containing  the  same  elements) . 

In  the  usual  way,  we  will  use  the  imbedding 
VC  H  =  H'  C  V'  , 


where  the  injections  from  V  into  H  and  from  H  into  V*  are  continuous  with 
dense  ranges.  We  denote  by  Ay  the  Riesz  map  from  V  onto  its  dual  V';  i.e.  , 


<v,vt>v  =  (A„v.)v,  v. , v  e  V. 

V  (2.5) 

Then  AQ  is  the  restriction  of  Ay  to  D(AQ)  in  the  sense  that 
(AyVjJv  =  <v,^)v1>H,  Vj^  e  D(A0),  v  e  V. 

(2.6) 


Now  we  define  the  total  energy  space  E  =  V  X  H,  noting  that  when  AQ  is 

•  2 

coercive  and  x(t)  is  the  solution  to  (2.1),  then  | | (x(t) ,x(t)) | |g  is  twice  the 
total  energy  (kinetic  plus  potential)  in  the  system.  We  want  to  write  (2.1) 
as  a  first-order  evolution  equation  on  E.  To  to  this,  we  must  determine  the 
appropriate  semigroup  generator  for  the  open-loop  system.  We  will  derive  this 
generator  by  constructing  its  inverse  explicitly,  and  then  we  will  try  to  con¬ 
vince  the  reader  that  we  do  have  the  appropriate  open-loop  semigroup  genera¬ 
tor.  The  approach  seems  mathematically  efficient,  and  we  will  need  the 
inverse  of  the  generator  for  the  approximation  scheme.  First,  we  must  state 


our  precise  hypotheses  on  damping  and  discuss  its  representation. 


Xhs  Hanging  Fang-tig  nal  and  Operator 

Actually,  we  do  not  require  an  operator  DQ  defined  from  some  subset  of  H 
into  H.  Rather,  we  assume  only  that  there  exists  a  damping  functional 


d0(vl’v2):  v  *  V-*R 


(2.7) 


such  that  dQ  is  bilinear,  symmetric,  continuous  (on  V  %  V)  and  nonnegative. 


If  we  have  a  symmetric,  nonnegative  damping  operator  DQ  defined  on  D(A0) 
such  that  DQ  is  bounded  relative  to  Aq,  then  ^qv^.vj^  defines  a  bilinear, 
symmetric,  bounded,  nonnegative  functional  on  a  dense  subset  of  V  X  V.  In 
this  case,  the  unique  extension  of  this  functional  to  V  x  V  is  dQ.  (That  DQ 
being  AQ-bounded  implies  continuity  of  with  respect  to  the  V-norm 

follows  from  [Kl,  Theorem  4.12,  p.  292].) 


Under  our  hypotheses  on  d0,  there  is  a  unique  bounded  linear  operator  AD 
from  V  into  V  such  that 


Vv.v.)  =  (Anv1)v,  v..v  c  V. 

U  1  (2.8) 

The  operator  (Ay^Ap)  is  then  a  bounded  linear  operator  from  V  to  V,  and 

(Ay^Ajj)  is  selfadjoint  (on  V)  because  dQ  is  symmetric.  Also 

d0(v,vi)  =  <v,Av1Adv1>v  =  <Ay1ADv,v1>y  v1>V  e  V. 

(2.9) 

Remark  2.2.  We  chose  to  begin  our  description  of  the  control  system  model 
with  (2.1)  because  its  form  is  familiar  in  the  context  of  flexible  structures. 
The  stiffness  operator  Aq,  for  example,  is  the  infinite  dimensional  analogue 
of  the  stiffness  matrix  in  finite  dimensional  structural  analysis.  In  appli¬ 
cations  like  the  example  in  Section  6,  though,  it  is  often  easier  to  begin 
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with  a  strain-energy  functional  from  which  the  correct  strain-energy  inner 
product  for  V  is  obvious.  The  stiffness  operator  is  defined  then  in  terms  of 
the  Riesz  map  for  V  (see  [S3]  for  this  approach),  rather  than  V  being  defined 


in  terms  of  the  stiffness  operator;  specifically,  AQ  is  defined  by  (2.6)  with 
**  -1 

D(Ag)  =  Ay  H.  Either  way,  the  relationship  between  Aq  and  V  is  the  same.  But 
the  only  thing  that  needs  to  be  computed  in  applications  is  the  V-inner  prcr- 
duct;  an  explicit  Aq  need  not  be  written  down.  O 


Ike  Semigroup  Generator 

We  define  A-1  e  L(E, E)  by 


A"1 


* -1a  <-l 

“Ay  Ajj  -Aq 

I  o 


(2.10) 


This  operator  is  clearly  one-to-one,  and  its  range  is  dense,  since  V  is  dense 
in  H  and  D(AQ)  ia  dense  in  V.  Now,  we  take 


Direct  calculation  of  the  inner  product  shows 


(2.11) 


(X-1 


v 

h 


vh>E 


=  -d0(h,h). 


(2.12) 


so  that  A  is  dissipative  with  dense  domain.  Also,  since  D(A  )  =  E,  A  is  max¬ 
imal  dissipative  by  [Gl,  Theorem  2.1].  Therefore,  t  generates  a  CQ- 
contraction  semigroup  on  E. 


Finally,  the  open-loop  semigroup  generator  is 


« 


« 
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A  + 


D(A) 


D  (A) 


where  A^  is  the  bounded  linear  operator  discussed  above.  With 
B  =  J  e  L(r“.E). 

the  first-order  form  of  (2.1)  is 

z(t)  -  Az(t)  +  Bu(t) ,  t>0 


where  z  =  (x. x)  e  e. 


(2.13) 


(2.14) 


(2.15) 


To  see  that  A  is  indeed  the  appropriate  open-loop  semigroup  generator, 
suppose  that  Aq  is  coercive  (so  that  A^  =  0)  and  that  we  have  a  symmetric, 
nonnegative  AQ-bounded  damping  operator  DQ.  Then  the  appropriate  generator 
should  be  a  maximal  dissipative  extension  of  the  operator 


o 

A 


o 

D(A)  =  D(A0)XD(Aq). 


(2.16) 


o 

It  is  shown  in  [Gl,  Section  2]  that  A  has  a  unique  maximal  dissipative  exten¬ 
sion,  and  it  can  be  shown  easily  that  the  A  defined  above  is  an  extension  of 
after  noting  that,  in  the  present  case. 


<avad)id(  -  aJJ  Dq. 

(2.17) 


We  should  note  that  Showalter  [S3,  Chapter  VI]  elegantly  derives  a 
semigroup  generator  for  a  class  of  second-order  systems  that  includes  the 
flexible- structure  model  here.  The  presentation  here  is  most  useful  for  our 
approximation  theory  because  of  the  explicit  construction  of  the  inverse  of 
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>o 


the  semigroup  generator.  For  the  purposes  of  this  paper,  we  do  not  need  to 
characterize  the  operator  A  itself  more  explicitly,  but  we  should  make  the 
following  points. 


First, 


from  A 


we  see 


D(A)={(x,x) :  x  e  V, 


_i 

x  +  Av  Adx  e  D(Aq)  )  . 


(2.18) 


In  many  applications,  especially  those  involving  beams,  the  "natural  boundary 
conditions"  can  be  determined  from  (2.18)  and  the  boundary  conditions  included 
in  the  definition  of  DAQ) .  In  the  case  of  a  damping  operator  that  is  bounded 
relative  to  AQ  ,  D(A)  =  D(Aq)xV.  If  the  damping  operator  is  bounded  rela¬ 
tive  to  Aq  for  p  <  1,  then  A  has  compact  resolvent. 


In  many  structural  applications,  the  open- loop  semigroup  is  analytic, 
although  this  has  been  proved  only  for  certain  important  cases.  Showalter 
obtains  an  analytic  semigroup  when  the  damping  functional  is  V-coercive;  for 
example,  when  there  exists  a  damping  operator  DQ  that  is  both  Astounded  and 
as  strong  as  AQ.  Such  a  damping  operator  results  from  the  Kelvin- Voigt 
viscoelastic  material  model.  Also,  it  can  be  shown  that  the  semigroup  is  ana¬ 
lytic  for  a  damping  operator  equal  to  c^Aq  for  1/2  i  |i  i  1  and  cQ  a  positive 
scalar.  The  case  p  =  1/2,  which  produces  the  same  damping  ratio  in  all  modes, 
is  especially  common  in  structural  models,  and  Chen  and  Russell  [Cl]  have 

shown  that  the  semigroup  is  analytic  for  a  more  general  class  of  damping 

1/2 

operators  involving  A q  . 


Finally,  we  can  guarantee  that  the  open-loop  semigroup  generator  is  a 
spectral  operator  (i.e. ,  its  eigenvectors  are  complete  in  E)  only  for  a  damp¬ 
ing  operator  that  is  a  linear  combination  of  an  H-bounded  operator  and  a 
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fractional  power  of  A^.  However,  nowhere  do  we  use  or  assume  anything  about 
the  eigenvectors  of  either  the  open-loop  or  the  closed-loop  semigroup  genera¬ 
tor.  The  natural  modes  —  of  undamped  free  vibration  —  in  (2.2)  are  always 
complete  in  both  H  and  V. 


2.2  The  Adjoint  of  A 


(A  X) 


Since  (Ay  AD)  is  selfadjoint  on  V,  direct  calculation  shows  that  A 

*  1 

—  the  adjoint  of  A  with  respect  to  the  E- inner  product  —  is 


T** 


A-1A  »f-l 

'Ay  Ap  Aq 

-I  0 


(2.19) 


Then  If*  =  0 [  *  )  1.  Having^  *  explicitly  facilitates  proving  strong  conver¬ 

gence  for  approximating  adjoint  semigroups. 


2.3  Exponential  Stability 


The  following  theorem  says  that,  if  there  are  no  rigid-body  modes  and  if 
the  damping  is  ooercive  (basically,  all  structural  components  have  positive 
damping),  then  the  open-loop  system  is  uniformly  exponentially  stable.  That 
the  decay  rate  given  depends  only  on  the  lower  bound  for  the  stiffness  opera¬ 
tor  and  the  upper  and  lower  bounds  for  the  damping  functional  is  essential  for 
convergence  results  for  the  approximating  optimal  control  problems  of  subse¬ 
quent  sections.  The  theorem  is  a  generalization  of  Theorem  6.1  of  [Gl]  to 
allcw  more  general  damping,  but  the  proof  is  entirely  different  and  much 
nicer.  The  current  proof  uses  an  explicit  Lyapunov  functional  for  the  homo¬ 
geneous  part  of  the  system  in  (2.15).  Recall  that  T(')  is  the  open- loop  semi¬ 
group,  with  generator  A,  and  E  is  the  total  energy  space  V  x  H. 
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Theorem  2  .|> .  Suppose  that  AQ  and  dQ  are  H-eoercive.  Let  p  be  the  positive 
constant  in  (2.4),  and  let  6Q  and  be  positive  constants  such  that 


60iivii5i  Vv'v)  v  e  v* 


(2.20) 


Then 


I  lT(t )  |  |  i  (1  +*==  +  6160/2)1/2  exp(-t/(^~  +^|.+  Sj) )  .t^O 
■  ^  0  *  , . 

Proof.  For  y  >  max{  ^  .  ^-) ,  define  Q  e  L(E)  as 


(2.21) 


Q  = 


(yl+Ay^p)  A^1 

I  Yl 


(2.22) 


Since  (AyMp)  is  selfadjoint  and  nonnegative  on  V,  Q  is  selfadjoint  and  coer¬ 
cive  on  E.  Define  the  functional  p(*)  on  E  by 


p(z)  =  <Q2.z>e  -  T ( I lx| |y  +  | |x| ly)  +  2<x,  x>h  +  dQ(x,x) , 


where  z  =  (x.x).  From  (2.4).  we  have 


(2.23) 


2 | <x, x>Hl  <  ^  l|x||y||x||H  i 


(2.24) 


so  that 


(rypJnzi|E  p(z)  1  p_1hzMe 


with 


0  =  (Y+  fT+  6i} ~1’ 


(2.25) 


(2.26) 
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Next  take  z  =  (x,x)  e  D(A)  and  set  A(x,x)  =  (y,y)  e  E,  or 


(x.x)  =  A-1(y»y)  =  (-Ay1ADy-A^1y,y> 


(2.27) 


Note  x  =  y  e  V.  Now, 


<QAz,z>e  =  <Q(y,y) ,A  1(y,y)>E 


-  Y<y*Ay1ADy>v  -  I  |Av_1ADy+A^1y|  |y  +  I  |y  1 1| 


=  -  l||x| |v  +  yd0(x.x) 


I  I*X|  1^3 


(2.28) 


From  (2.20)  then. 


<qaz,z>e  i  -  [||x|  |y  +  (yV1)MxIIh1  i  “I  I* I  l|« 


(2.29) 


Therefore  p(’)  is  a  Lyapunov  functional,  and  the  theorem  follows  from  (2.25), 
(2.26)  and  (2.29),  with  y  =  y^.+  .  □ 


•  • 


•  • 
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3.  The  Optimal  Control  Problem 

Subsection  3.1  presents  some  preliminary  definitions  and  results  for  the 
optimal  linear-quadratic  regulator  problem  on  an  arbitrary  real  Hilbert  space. 
These  results  are  generic  in  the  sense  that  the  Hilbert  space  E  is  not  neces¬ 
sarily  the  energy  space  of  Section  2,  and  the  operators  A,  B,  etc.,  do  not 
necessarily  represent  an  abstract  flexible  structure  as  in  Section  2.  In  the 
second  half  of  the  paper,  having  such  generic  results  will  allow  us  to  obtain 
the  approximation  theory  for  the  infinite  dimensional  state  estimator  from  the 
analogous  results  for  the  control  problem.  Subsection  3.2  gives  some  impor¬ 
tant  implications  of  the  general  results  for  the  case  where  the  control  system 
is  that  defined  in  Section  2. 

3.1  The  Generic  Optimal  Regulator  Problem 

Let  a  linear  operator  A  generate  a  ^-semigroup  T(t)  on  a  real  Hilbert 
space  E,  and  suppose  B  e  L(Rm.E),  Q  e  L(E,E)  and  R  e  L(Rm) ,  with  Q  nonnegative 
and  selfadjoint  and  R  positive  definite  and  symmetric.  The  optimal  control 
problem  j&q  E  is  to  choose  the  control  u  e  (0,<»;Rm)  to  minimize  the  cost 
functional 


J(z(0),u)  -  /  ( <Qz(t)  ,z(t)  >-,  +  <Ru(t)  ,u(t)>nm)dt, 

0  < 

where  the  state  z(t)  is  given  by 

t 

z ( t )  =  T(t)z{0)  +  /T(t-n)Bu(n)dn,  t  2  o. 

0  <; 

Definition  3.1.  A  function  u  e  L2(0,®;U)  is  an  admissible  control  for  the 


22 


initial  state  z,  or  simply  an  admissible  control  for  2t>  if  J(z,u)  is  finite; 
i.e.  ,  if  the  state  z(t)  corresponding  to  the  control  u(t)  and  the  initial  con¬ 
dition  z(0)  =  z  is  in  L2(0,<=;E) . 


Definition  3  .2  Let  the  operators  A,  B,  Q,  and  R  be  as  defined  above.  An 
operator  II  in  L(E)  is  a  solution  of  the  Riccati  algebraic  equation  if  II  maps 
the  domain  of  A  into  the  domain  of  A*  and  satisfies  the  Riccati  algebraic 
equation 


A*n  +  n  a  -  nBR1B*n  +  q  =  o. 


(3.3) 


Theorem  3  .3  (Theorems  4.6  and  4.11  of  I G4 ]  ) .  There  exists  a  nonnegative  sel- 
fadjoint  solution  of  the  Riccati  algebraic  equation  if  and  only  if,  for  each 
z  e  E,  there  is  an  admissible  control  for  the  initial  state  z.  If  II  is  the 
minimal  nonnegative  selfadjoint  solution  of  (3.3),  then  the  unique  control 
u(  )  which  minimizes  J(z,u)  and  the  corresponding  optimal  trajectory  z(  )  are 
given  by 


u(t)  =  -R  *B  Ilz(t) 


(3.4) 


and 


z(t)  =  S(t)z, 

where  S(‘)  is  the  semigroup  generated  by  A-BR  *B^I.  Also, 

J(z,u)  =  min  J(z,v)  =  <  II  z,z>g 

v 

If,  for  each  initial  state  and  admissible  control. 


(3.5) 


(3.6) 


1 In  |  |z(t)|  | 
t— 


0 


(3.7) 


there  exists  at  most  one  nonnegative  selfadjoint  solution  of  (3.3).  If  Q  is 
coercive,  (3.7)  holds  for  each  initial  state  and  admissible  control  and  S(  •  ) 
is  uniformly  exponentially  stable.  □ 


We  will  refer  to  T(-)  as  the  open-loop  semigroup  and  to  S(*)  as  the 
optimal  closed-loop  semigroup. 

To  prepare  for  the  convergence  analysis  in  Sections  5  and  9,  we  must 
present  now  some  rather  arcane  estimates  for  the  decay  rate  of  the  closed-loop 
system  in  the  optimal  control  problem. 


Theorem  3 .4  Suppose  that  the  open-loop  semigroup  T(*)  satisfies 

a.  t 

I  I T ( t )  |  |  1  Mje  X  .  t  2  0, 


(3.8) 


for  positive  constants  and  c^,  that  II  is  the  minimal  nonnegative  selfad¬ 
joint  solution  to  (3.3),  and  that  S(t)  is  the  optimal  closed-loop  semigroup  in 
Theorem  3.3.  If  there  exists  a  constant  such  that,  for  each  z  e  E, 

QO 

/  1 1 S  ( t )  z  |  f  2dt  1  Hj(<nz,z>E  +  ||z||2) 

0 

9 

and  a  constant  Mq  such  that 

ll  n  ii  <  v 

then  there  exist  positive  constants  Mj  and  02 ,  which  are  functions  of  Mq, 

Mj,  and  only,  such  that 


(3.9) 


(3.10) 
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1 ls<t ) | |  i  Mje  2  ,  t  1  0. 


(3.11) 


Proof.  This  follows  easily  from  Theorems  2.2  and  4.7  of  [  ].  □ 


T.ftinmfl  3.§_.  Suppose  that  there  exist  positive  constants  H  and  a  such  that 


1 1 T ( t )  |  |  i  Me  at,  t  1  0. 


(3.12) 


If  z(0)  e  E,  h  e  L-(0,<°;E),  and  z(t)  =  T(t)z(0)  +  /  T(t-s)h(s)ds,  then 
2  0 


/||z(t)||Z  dt  1 
0 


|  |z(0>  I  I  +  |h|  |, 


(3.13) 


Proof.  The  result  follows  from  (3.12),  the  convolution  theorem  [Dl,  page  95l] 
and  the  triangle  inequality.  O 

3.6.  Suppose  that  E  is  finite  dimensional  and  that  the  pair  (Q.A)  is 
observable  (in  the  usual  finite  dimensional  sense).  Then  there  exists  a  con¬ 
stant  M,  which  is  a  function  of  A,  B  and  Q  only,  such  that 


/ 1  |z(t )  1 1 2dt  L  M  (/  (<Qz(t),z(t)>s.  +  |  |u(t)|  |2)dt, 

0  0  E 


(3.14) 


where  z(t)  is  given  by  (3.2). 


Proof.  The  proof,  which  is  at  most  a  mild  challenge,  is  based  on  the  fact 
that  the  observability  grammian 
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t  ATf  A*. 

W(t)=  /  eA  tQeAtdt 
0 

is  coercive  for  any  positive  t.  □ 


(3.15) 


The  next  theorem  says,  among  other  things,  that  if  the  open-loop  control 
system  decouples  into  a  finite  dimensional  part  that  is  stabilizable  (in  the 
usual  finite  dimensional  sense)  and  an  infinite  dimensional  part  that  is  uni¬ 
formly  exponentially  stable,  then  the  entire  system  is  uniformly  exponentially 
stabilizable,  so  that  (3.3)  has  a  nonnegative  selfadjoint  solution. 


Theorem  3.2*  Suppose  that  there  exists  a  finite  dimensional  subspace  E^, 
D(A)  such  that  EQ  and  E^*" reduce  A  (and  T(t)),  and  write 


[An  o  I  n 

rBn  i  n 

r°ii 

Q12 

A  = 

l  0  hll  = 

l  B2ll ' 

CS 

*  H 

a 

_ f 

^22 

(3.16) 


where  A^  and  A22  are  the  restrictions  of  A  to  EQ  and  D(A)  X  E^,  respectively. 
Simll arly. 


T(t) 


i 


rTu(t)  o 

L  0  T22^ 


(3.17) 


Also,  suppose  that  the  pair  (A^.Bjj)  is  stabilizable  and  that  there  exist 

9  9 

positive  constants  Mj,  and  0  such  that 


9 

,  -a,  t 

llT22(t)||  <  ,  t  2  0 


and 


(3.18) 


I 
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max{| |B| | , | |Q| | }  1  p. 


(3.19) 


i)  Then  there  exists  F  e  L(E,Rm)  such  that  A-BF  generates  a  uniformly 
exponentially  stable  semigroup  on  E.  Also,  (3.3)  has  a  nonnegative  selfad- 
joint  solution,  and  the  minimal  such  solution  satisfies  (3.10)  with  a  func- 

9  9 

tion  of  A11#  B1;l,  R,  and  p  only. 


ii)  If  Q12  =  0  and  the  pair  (Qj^.  ^j)  is  observable,  then  there  exists  a 
unique  nonnegative  selfadjoint  solution  II  to  (3.3),  and  there  exist  positive 

9  9 

constants  M2  and  —  which  depend  on  A^,  B^,  0^,  R,  M^,  and  p  only  — 
such  that  the  optimal  closed-loop  semigroup  satisfies 


l|S(t)||  1  Mje  °2t,  t  2  0. 


(3.20) 


Proof .  1)  To  say  that  <A11.B11)  is  stabilizable  means  that  there  exists  a 

linear  operator  F^  from  ^  to  such  that  each  eigenvalue  of  A^j-B^jFj^  has 
negative  real  part.  Hence  A-BF  generates  a  uniformly  exponentially  stable 
semigroup  if  F  =  [F21  0],  so  that  there  exists  an  admissible  control  for  each 
initial  condition  in  E. 


It  is  easy  to  write  down  an  upper  bound  for  the  performance  index  in 

9  9 

(3.1)  in  terms  of  R,  p  ,  M^,  o2  and  the  decay  rate  of  exp(  [A^  _B11  F11  ]t)* 

V  9  9 

That  the  in  (3.10)  depends  only  on  A11#  B21,  R,  M^,  and  P  follows  then 
from  the  fact  that  F21  is  a  function  of  A^  and  B^j. 

9  9 

ii)  Clearly,  (3.8)  holds  with  M2  and  a2  depending  only  on  Ajj,  B11#  Mj,  <*2 
and  F12»  Therefore,  we  have  (3.8)  and  (3.10)  with  the  bounds  depending  only 

9  9 

on  A11#  Bj^j,  ,  R,  Mj ,  02  and  P  .  Finally,  the  existence  of  an  ^  for 
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(3.9)  which  depends  only  on  these  parameters  follows  from  using  (3.1)  and 
(3.6)  in  applying  Lemma  3.6  to  the  part  of  the  system  on  EQ  and  then  Lemma  3.5 
to  the  part  of  the  system  on  E^-.  Part  ii)  of  this  theorem  then  follows  from 
Theorem  3 .4 .  □ 

9 

Remark  £.£  When  we  say  in  Theorem  3.4  that  Mj  and  a2  are  functions  of  Mo* 

M^,  and  only,  we  mean,  for  example,  that  for  two  optimal  control  problems 
on  different  spaces  E,  with  different  operators  A,  B,  etc. ,  if  the  same  con- 

9 

stants  V  V  «!•  and  a1  work  in  (3.8)-(3.10)  for  both  problems,  then  the 
same  constants  Mj  and  will  work  in  (3.11)  for  both  problems.  Similarly,  in 

9  9 

Theorem  3.7  ii) ,  as  long  as  EQ,  A11<  B^,  R,  Mj,  and  p  remain  the 

same,  the  same  M2  and  02  will  work  in  (3.20)  even  if  E^*,  and  ^22 

change.  □ 

3.2  Application  to  Optimal  Control  of  Flexible  Structures 

For  the  rest  of  this  section,  A^,  Aj,  A,  T(t),  Bq  and  B  are  the  operators 
defined  in  Section  2.1,  and  E  =  V  X  H  is  the  energy  space  defined  there. 

Rama r lr  3.9.  Theorem  3.7  is  useful  mainly  when  all  but  a  finite  number  of 
modes  have  coercive  damping  in  the  open- loop  system  and  the  damped  and 
undamped  parts  of  the  open- loop  system  remain  orthogonal.  This  is  the  case, 
for  example,  with  modal  damping.  The  next  theorem  does  not  require  ortho¬ 
gonality  of  the  damped  and  undamped  parts  of  the  system,  but  it  does  require 
an  independent  actuator  for  each  undamped  mode.  The  situation  of  Theorem  3.10 
is  typical  in  aerospace  structures:  Any  elastic  component  should  have  some 
structural  damping,  but  rigid-body  modes  are  common;  for  a  structure  to  be 
controllable,  an  actuator  is  required  for  each  rigid-body  mode.  □ 


XfafiflESB  1.10.  i)  Suppose  that  =  bqbJ  and  that  TQ  =  Aq  +  and  dQ  =  dQ  + 
are  H-coercive,  so  that  there  exist  positive  constants  p  ,  y  and  p  such 
that,  for  all  v  e  V, 


llvllj  1  PllvjlJ, 


dn(v»v)  2  p||v||„. 


d0(v,v)  i.  rlMly, 


(3.21) 


(3.22) 


(3.23) 


max{llB0||,||Q||.||R||}  i  p 


(3.24) 


(The  V-continuity  of  dQ  implies  (3.23).)  Then  (3.3)  has  a  minimal  nonnegative 

9 

selfadjoint  solution  II.  which  satisfies  (3.10)  with  a  function  of  p,  y 
and  p  only  . 

ii)  Suppose  also  that 


<0z,z>E  2  p| I z | |g,  z  e  E. 


Then  the  optimal  closed-loop  semigroup  satisfies 


IIS(t)||  i  ^e  ,  t  2  0, 


where  and  a-  are  positive  constants  depending  on  p,  y  and  p  only. 


Proof,  i)  The  suboptimal  control 


(3.25) 


(3.26) 
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(3.27) 


u(t)  =  -BQ[x(t)  +  x(  t ) ] 

produces  a  closed-loop  system  with  exponential  decay  at  least  as  fast  as  that 
in  Theofem  2.1.  The  required  upper  bound  in  (3.10)  follows  then  from  (3.1), 
(3.6)  and  (3.24) 

ii)  In  this  case,  the  Mp  in  (3.9)  is  p,  and  Theorem  3.4  yields  the  result.  □ 


Now  we  will  consider  the  structure  of  the  optimal  control  law  In  more 
detail.  Since  II  e  L(E,E)  and  E  =  V  X  H,  we  can  write 


n  = 


Do  n1 
nj  n. 


(3.28) 


where  Dp  e  L  ( V,  V)  ,  s  T  H,V)  ,  s  L(H,  H),  and  Dp  and  II,  are  nonnegative 
and  selfadjoint.  With  z  =  (x.x),  as  in  Section  2,  (3.4)  becomes 


u(t)  =  -R_1Bp [  D^x(t )  +  D^U)] 


(3.29) 


Since  Bp  e  L(Rm,H),  we  must  have  vectors  b^e  H,  1£  i  1  m,  such  that 


a 


Bou  =  ,z.biui 


i=l 


(3.30) 


for 


U  =  [u.u-  u  ]Te  Rm. 

±2  nj 

Aiso,  for  h  e  H, 

Boh  =  [<ba,h>H  <b2,h>H  *•*  <bm,h>H]T. 

* 

Since  ILx(t)  and  ILx(t)  are  elements  of  H,  we  see  from  (3.29)  and 


(3.32) 


(3.32)  that  the  components  of  the  optimal  control  have  the  feedback  form 


where  fA 


Uj_(t)  =  -  <fi,x(t)>v  -  <gi,x(t)>H,  i  =  l,...,m, 

e  V  and  e  H  are  given  by 


(3.33) 


®  -1 

i  (r  n 

j=i 


U 


nlbj. 


(3.34a) 


g<  "  1  (R  1)ij  ^2bj, 


1=1. 


i  m. 


j=l 


(3.34b) 


We  call  f  and  g 


gains. 


4.  The  Approximation  Scheme 

4.1  Approximation  of  the  Open- loop  System 

Hypothesis  ±.1.  There  exists  a  sequence  of  finite  dimensional  subspaces  Vn  0f 
V  such  that  the  sequence  of  orthogonal  projections  Pyn  converges  V-strongly  to 
the  identity,  where  Pyn  is  the  V-projection  onto  Vn.  Also,  each  Vn  is  the 
span  of  n  linearly  independent  vectors  e^. 

Since  it  should  cause  no  confusion,  we  will  omit  the  subscript  n  and  write 
just  keeping  in  mind  that  the  basis  vectors  may  change  from  one  Vn  to 

another,  as  in  most  finite  element  schemes.  Also,  we  will  refer  to  the  Hil¬ 
bert  space  En  =  Vn  X  Vn>  which  has  the  same  inner  product  as  E  =  V  X  H. 

For  n  2  1,  we  approximate  x(t)  by 

n 

xn(t)  *  I  5 .i(t)e,, 

J=1  -  J 

where  S(t)  *  (^(t).  52(t),  ••• 

MnV<t)  +  Dn?(t)  + 

and  the  mass  matrix  Mn,  damping  matrix  Dn,  stiffness  matrix  Kn  ,  and  actuator 
influence  matrix  are  given  by 


(4.1) 

$n(t))T  satisfies 
K^Ct)  =  B?u(t), 

(4.2) 
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M°  =  [<e1,eJ>H],  Dn  =  [d^e^)], 

=  I<ei*ej>v1  ~  [<A1ei,ej>Hl, 

B0  =  * 

Of  course,  (4.2)  can  be  written  as 
n  =  Ann  +  Bnu 

where 

H  =  UT.  «T ]T 

and 


(4.3) 


(4.2*) 


(4.4) 


Note:  Throughout  this  paper,  we  use  the  superscript  n  in  the  designation  of 
matrices  in  the  ntb  approximating  system  and  control  problem,  like 
An,  Bn,  M0,  etc.  Hence  the  superscript  n  indicates  the  order  of 
approximation  —  and  it  never  indicates  a  power  of  the  matrix.  By  M-n, 
we  denote  the  inverse  of  the  mass  matrix  Mn. 

In  the  designation  of  a  linear  operator  in  the  nth  approximation,  we  use  the 
the  subscript  n.  For  example,  Afl  and  BR  are  the  operators  whose  matrix 
representations  are  An  and  Bn,  respectively. 
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For  convergence  analysis,  it  is  useful  to  note  that  (4.1)  and  (4.2)  or 


(4.2')  are  equivalent  to 


zn(t)  =  Vn(t)  +  Bnu(t)' 

(4 

where  z  =  (x  .x  )e  E  .  and  e  L(E  )  and  e  L(Rm.  E  )  are  the  operators 
n  n  n  n  nn  n  n 

whose  matrix  representations  are  given  in  (4.5).  Also,  for  any  real  A, 


is  equivalent  to 

(AV1  +  ADn  +  Kn)a1  =  (AMn  +  D11)/)1  +  M1^2 

and 

a2  =  Aa1  -  p1 


(4.7) 


(4.8) 


(4.9) 


vJ  = 
n 

(Substituting  An 


Z  a^e  and  h^  =  I  p!?e, ,  J  =  1»2. 

i=l  1  n  i=i  A 

and  (4.10)  into  (4.7)  yields  (4.8)  and  (4.9)). 


(4.10) 


Next,  we  will  prepare  to  invoke  the  Trotter-Kato  semigroup  approximation 
theorem  to  show  how  (4.2),  (4.2')  and  (4.6)  approximate  (2.1)  and  (2.15). 
First,  we  will  treat  the  case  in  which  Aq  is  coercive  (no  rigid— body  modes), 
so  that  A^=0  and  Aq=Aq;  the  general  case  Is  a  straight-forward  extension.  For 
Aq  coercive,  the  open-loop  semigroup  generator  A  is  maximal  dissipative.  Also, 
for  each  n,  A^  is  dissipative  on  En<  The  main  idea  here  is  to  project  (A-A)  1 
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onto  Vn  in  a  certain  inner  product  and  observe  that  the  result  is  exactly 
(X-Afi)  *,  where  Ar  is  the  operator  on  in  (4.6)  and  (4.7).  Of  course,  we 
need  only  do  this  for  real  X>0. 


For  real  X>0  then,  define  an  innerproduct  on  V  by 


<*.*>,  =  X2<*,*>„  +  Xdn(\*>  + 

K  M  °  V  (4.11) 

Under  the  hypotheses  in  Section  2  on  dQ,  <‘»‘\  is  norm- equivalent  to  <*,*>y. 

For  n  2  1,  let  pn(X)  be  the  projection  of  V  onto  VR  in  the  inner  product 


Now  let  bjjhj  e  H  and  note  that 


(X-A) 


is  equivalent  to 


£]*(S 


(4.12) 


(»2)‘  rl[^)‘(2»)l  • 

With  A-1  from  (2.10),  (4.13)  is  equivalent  to 


(4.13) 


(I  +  XA^Ap  +  X2Aq1)v1  =  (XA01  +  Av1AD)h1  +  A()1h2 


(4.14) 


and 


v2  =  Xv1  -  h1 . 

If 

vj  -  and  v2  .  PnU),2. 

it  follows  from  (4.11)  and  (4.14)  that 


(4.15) 


(4.16) 
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<ei*  vn>X 


<ei'vl>X 


~  ^2<ei»A0A^lvl>H  +  X<ei,Ay1ADv1>v  +  <ei,v1>^ 


(4.17) 


=  <e1,(X2A^1+XA“1AD+I)v1>v 


=  <ei(XAo1+Ay1AD)h1+AQ1h2>v. 


and  from  (4.15)  that 


<e.,v„\ 
i  n  x 


<ei.v2>x 


X<ei,v  >x 


<ei*h  >X 


(4.18) 


Now,  for  h1  =  h2  s  VR,  h2  *  h2  e  Vn>  and  v2,  v2,  h2  and  h2  written  as  in 
(4.10),  (4.17)  and  (4.18)  yield  (4.8)  and  (4.9)  again. 


This  shows  that 


'P_(X)  n  1  -1  -1 

0  Pn(X)|U~A)  *En  =  (X_V  ' 


which  yields 


[Pn(X)  0  1  -1 

[  0  Pn(X)  ]U~A) 


PEn  =  U_An)  PEn 


(4.19) 


(4.20) 


where  PF  the  E~1 


of  E  onto  En.  The  projection  PEn  can  be  written 


[PVn  0  I 
=  l  0  PHnl  ' 


(4.21) 


where  PVn  is  the  V- projection  onto  Vn>  as  before,  and  PHn  is  the  H-projection 
onto  Vn>  Since  the  V-norm  is  stronger  than  the  H-norm  and  the  norm  induced,  by 
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the  X-inner product  is  equivalent  to  the  V-norm,  it  follows  from  Hypothesis  4.1 
that  (X-An)  1  PEn  converges  E- strongly  to  (X-A)  1  as  n-»  ®.  Now,  with  A^ 
extended  to  Ej^  as,  say,  n(PEn-I) ,  Trotter-Kato  [Kl,  page  504,  Theorem  2.16] 
yields  the  following. 

Theorem  4.2.  For  Aq  coercive,  let  Tn  (')  be  the  (contraction)  semigroup  gen¬ 
erated  on  En  by  AR.  Then,  for  each  t  10,  ^n(t)Pg  °°nver6es  strongly  to 
T(t),  uniformly  in  t  for  t  in  bounded  intervals. 

In  the  general  case,  when  A^  is  not  coercive,  the  open-loop  generator  A 
is  obtained  from  the  dissipative  *k  by  the  bounded  perturbation  in  (2.13),  30 
that  [G3,  Theorem  6.6]  yields  the  following  generalization  of  Theorem  4.2. 

£oro,Uary  4*1 .  Let  Tn(’)  be  the  semigroup  generated  on  E^  by  Afi.  Then,  for 
each  t  2  0,  TQ(t  )PEn  converges  strongly  to  T(t),  uniformly  in  t  for  t  in 
bounded  intervals. 

Theorem  ±.±.  When  A  has  compact  resolvent,  (X-An)-1PEn  converges  in  L(E)  to 
(X-A)-1. 

Proof .  This  follows  from  (4.20)  and  a  standard  result  that  the  projections  of 

a  compact  linear  operator  onto  a  sequence  of  subspaces  converge  in  norm  if  the 

projections  converge  strongly  to  the  identity,  as  do  P^  and  P  (X).  □ 

bn  n 

That  the  adjoint  semigroups  also  converge  strongly  follows  from  an  argu¬ 
ment  entirely  analogous  to  the  proof  of  Theorem  4.2.  In  particular,  equations 
like  (4.11)-(4.17)  are  used  to  show  that 
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P  (A)]  a_A  }  PEn  “  (X_An>  lpEn 
n  J 


(4.22) 


In  showing  this,  A  *  is  used  as  A-*  was  used  above.  Also,  care  oust  be  taken 

* 

to  calculate  Ar  with  respect  to  the  E-inner  product.  The  result  is 

Theorem  4 .5 .  Let  T^C)  be  the  sequence  of  semigroups  in  Corollary  4.3.  Then, 
for  each  t  2  0,  T  (t)P.,  converges  strongly  to  T  (t),  uniformly  in  t  for  t  in 

n  tin 

bounded  intervals. 


Finally,  for  the  approximation  to  the  actuator  influence  operator  B  e 
L(Rm,E),  recall  Bn  e  L(Rm»  En) ,  the  operator  whose  matrix  representation  is 
the  matrix  Bn  in  (4.5).  From  (4.3),  it  follows  that 


Bn  PEnB* 


(4.23) 


Since  B  has  finite  rank  m,  Bn  and  B*  converge  in  norm  to  B  and  B*,  respec¬ 
tively. 
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4.2  The  Approximating  Optimal  Control  Problems 

x.u  • 

The  ntn  optimal  control  problem  is:  given  zn(0)  =  (xn(0) ,xn(0) )  e  En. 
choose  u  e  L^O.-iR0)  to  minimize 


00 

Jn*2n^'u^  =  H  <Qnzn(t)  »zn(t)>g  +  <Ru(t)  ,u(t)  >Rm)dt, 

0  (4.24) 

where  Qn  =  PgnQ|gn»  We  assume: 


Hypothesis  4.6.  For  each  n  2  1  and  zn(0)e  En,  there  exists  an  admissible  con¬ 
trol  (Definition  3.1)  for  (4.6)  and  (4.24) .  □ 


A  sufficient  condition  for  Hypothesis  4.6  is  that,  for  each  n,  the  system 
(An,Bn)  be  stabilizable. 


By  Theorem  3.1,  the  optimal  control  uR(t)  has  the  feedback  form 


u„(t>  - 


(4.25) 


where  is  a  linear  operator  on  En,  Dn  is  nonnegative  and  self ad joint,  and 
satisfies  the  Riccati  equation 


AX  *  V.  -  W*»X  *  °n  ■  »• 

( 4  «2* 

A3  a  result  of  Hypothesis  4.6,  (4.26)  has  at  least  one  nonnegative,  selfad- 
joint  solution.  The  minimal  such  solution  is  the  correct  here.  If  the 
system  (An,Qn)  is  observable,  then  1^  is  the  unique  nonnegative,  selfadjoint 
solution  to  (4.26),  and  is  positive  definite.  If  we  write  as 
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tj  ^On  Din 

^  =  nln 


(4.27) 


then  (4.25)  becomes 


vt)  =  - 


lOtninxn(t)  +  ^n*n(t)] 


(4.28) 


The  feedback  law  (4.28)  can  be  written  in  functional-feedback  form,  just 
as  in  Section  3.  We  have 


Un(t)  =  [uln(t>  u2n(t)-*-umn(t)1 


(4.29) 


where 


uin(t)  =  '  <fin’xn(t)>V  '  <«in'xn(t)>H'  1  <  i  <  “• 


(4.30) 


I  (R  ^  DinPHnb.j, 


1  1  i  i  m. 


)ijI^nPHnbj  '  1  1  i  <  m. 


(4.31a) 


(4.31b) 


Of  course,  f  and  gin  are  the  ntb  approximations  to  the  functional  gains  f^ 
and  g±  in  (3.34)  . 

In  Section  5,  (4.25)-(4.31)  will  be  useful  for  studying  how  the  solution 
to  the  nbb  optimal  control  problem  converges  to  the  solution  to  the  original 
problem  of  Section  3,  but  for  numerical  solution  of  the  nth  problem,  we  need 
the  matrix  representations  of  these  equations. 

We  will  need  the  following  grammian  matrices: 
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1 


(4.32) 


wn  =  P”  «]. 

LO  MnJ 


(4.33) 


Note:  The  matrix  W-n  will  be  the  inverse  of  Wn.  The  superscript  n  on  any  matrix 
indicates  the  order  of  approximation,  not  a  power  of  the  matrix.  Also,  recall 
the  note  following  (4.5) . 

Now  recall  Qn  =  PEnQIEn.  Since  Q  =  Q*  e  L(E)  and  E  =  V  JC  H,  we  can  write 


*0 

°2 


where  QQ  =  C^e  L(V),  e  L(H,V),  and  Q2  = 
lation  shows  that 


Qn  =  w-nQn’ 


(4.34) 


Q2  e  LOO  .  Straightforward  calcu- 


(4.35) 


where  Qn  is  the  matrix  representation  of  Qn  and  Qn  is  the  nonnegative,  sym¬ 
metric  matrix 


oj  qJ 

of  Qj 


(4.36) 


Q0  ^<ei’C!Oej>V-*  ' 


t ^6^ , Q^e j )y]  , 


(4.37) 


°2  [<ei»02ej>H] * 

Also,  recall  that  A^  and  Bn  are  the  operators  whose  matrix  representations  are 
given  by  (4.5),  and  note  that  the  matrix  representations  of  A*  and  B*  are 
W-n(An)TWn  and  (Bn)TWn,  respectively. 


With  the  matrix  representation  of  1^  denoted  by  if#  the  Riccati  operator 
equation  (4.26)  is  equivalent  to  the  Riccati  matrix  equation 


W~n(An)TWn  If1  +  IPA11  -DnBnR-1(Bn)TWnIIn  +  Qn  =  0. 

(4.38) 

While 

is  self  ad  joint,  if  in  general  is  not  symmetric,  but  the  matrix 


fr1  =  WHlf1 

(4.39) 

is  symmetric  and  nonnegative,  and  positive  definite  if  Hn  is.  Premultiplying 

(4.39)  by  Wn,  we  obtain 


(An)T  if  +  if  An  -  fr1BnR_1(Bn)Tnn  +  Qn  =  0, 

which  is  the  Riccati  matrix  equation  to  be  solved  numerically. 


(4.40) 


Now  we  need  one  more  set  of  matrix  equations  for  the  numerical  solution 
of  the  nth  optimal  control  problem.  Since  the  functional  gains  f^  and  gin 
are  elonents  of  V  ,  they  can  be  written  as 


n 

f.  =  7  B.  e. 

in  jrj  po  0 

and 

«in  = 

n  gi 

Z  ^j  ej’ 
j=l  J  J 

i  =  1, . . 

.  ,m, 

(4.41) 

fi  gi  „ 

where  3  ,  p  e  R.  We  reed 

equations 

fi 

for  p 

gi 

and  p  in 

terms  of  If1. 

One 
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way  to  get  these  equations  is  to  partition  IP  (obtained  from  (4.39))  and  then 
work  out  the  matrix  representation  of  (4.31).  However,  another  approach  is 
more  instructive  because  it  relates  the  present  Hilbert  space  methods  to  the 
standard  finite  dimensional  solution  of  the  n^  optimal  control  problem. 


The  n^h  optimal  control  problem  can  be  stated  equivalently  as:  given 
ij(0)  =  [£(0)T,£(0)T]  e  R2n  choose  u  e  L2(0,»jRm)  to  minimize 

00 

J  (tj(0)  ,u)  =  J  (T\(t)T5rxn(t)  +  u(t)TRu(t))dt, 

11  0  (4 

where  n(t)  =  [£(t)T,£(t)^lT  satisfies  (4.2’).  For  (4.2')  and  (4.42),  the 

optimal  control  law  is 


U  =  -R-VtfVt) 


(4.43) 


where  IP  is  the  minimal  nonnegative,  symmetric  solution  to  (4.40). 


Since  5  is  related  to  xr  by  (4.1),  the  optimal  control  un  in  (4.43)  must 
be  equal  to  the  optimal  control  ufi  in  (4.29)-(4.31) .  Substituting  (4.41)  into 
(4.30)  yields 


uin  =  -  (pfi)TKn5  -  (pgi)TMn£  =  -C(pV  (pgl)T]wnn. 


(4.44) 


Then,  using  (4.44)  and  equating  (4.29)  and  (4.43)  yields 


fl  f2  ... 
P  P 

8l  go 

LP  p 


‘m 


P 


=  w_nnnBnR  1 


(4.45) 


W§  now  have  the  complete  solution  le  the  nth  optimal  control  problem: 


4 
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The  Riccati  matrix  equation  (4.40)  is  solved  for  If1;  then  the  optimal  control 
is  given  by  (4.43),  and  equivalently  by  (4.29)-(4.30)  with  the  functional 
gains  f^n  and  g^n  given  by  (4.41)  and  (4.45).  In  the  next  section,  we  will 
give  sufficient  conditions  for  the  solution  to  the  nfch  optimal  oontrol  problem 
to  oonverge  to  the  solution  to  the  optimal  control  problem  in  Section  3  for 
the  original  infinite  dimensional  system. 
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5.  Convergence 


As  in  Section  3,  subsection  5.1  will  state  some  results  for  the  optimal 
linear  regulator  problem  involving  generic  linear  operators  A,  B,  Q,  etc. ,  o 
an  arbitrary  real  Hilbert  space  E,  and  subsection  5.2  will  expand  upon  these 
results  for  the  particular  class  of  control  problems  treated  in  this  paper. 


5.1  Generic  Approximation  Results 


Let  the  Hilbert  space  E  and  the  linear  operators  A,  T(*),  B,  Q  and  R  be 

as  in  Section  3.  Suppose  that  there  is  a  sequence  of  finite  dimensional  sub 

spaces  En,  with  the  projection  of  E  onto  En  denoted  by  Pgn,  such  that  Pgn  co 

verges  strongly  to  the  identity  as  n  4  «  ,  and  suppose  that  there  exist 

sequences  of  operators  A  8  L(E  ) ,  B„  b  L(Ra,E  ),  Q  =  Q*  e  L(E  ),  <3  2  0 

n  nn  nnn  nn 

such  that  we  have  the  following  strong  convergence.  For  all  z  e  E  and  t.^0. 


exp(Ant)PEnz  -4  T(t)z 

and 


(5. 


exp(A*t)PEnz  -4  T*(t)z 

as  n  4  •,  uniformly  in  t  for  t  in  bounded  intervals;  for  each  u  e  Rm, 
Bnu  -4  Bu; 
for  each  z  e  E, 


(5.: 


(5.: 


QnPEns  Qz> 


(5.- 


Theorem  5.1.  Suppose  that  for  each  n  there  is  a  nonnegative,  selfadjolnt 
linear  operator  IL  on  En  which  satisfies  the  Riccati  algebraic  equation 


‘X  *  V.  -  W1bX  *  °n  -  »•  (J  , 

If  there  exist  positive  constants  M  and  0,  Independent  of  n,  such  that 

||exp([An-BnR' Vnn]t)||  <  Me"pt.  t  2  0. 

(5.6 

and  if  ||  1 1  is  bounded  uniformly  in  n,  then  the  Riccati  algebraic  equation 

(3.3)  has  a  nonnegative  selfadjoint  solution  II  ,  and,  for  each  z  e  E, 


UnW  132 

and 


(5.7) 


e*»('*„-y"lBXu>pEn‘  ->  S(t)* 

(5*0/ 

uniformly  in  t  2  0,  where  SC)  is  the  semigroup  generated  by  A-BR‘  _1b1i  .  If 
there  exists  a  positive  constant  6,  independent  of  n,  such  that 

Qn  ^  8' 

(5.9) 

then  ||iy,  being  bounded  uniformly  in  n  guarantees  the  existence  of  positive 
constants  M  and  0  for  which  (5.6)  holds  for  all  n. 


Proof.  The  theorem  follows  from  Theorem  5.3  of  [G4]  when  the  operators  An>  Qn 
and  Hjj  are  extended  to  all  of  E  by  defining  them  appropriately  on  EjJ-.  For  the 
details  of  this  procedure,  see  Section  4  of  [Gl] .  Or  better.  Banks  and  Kun- 
isch  [B6]  have  modified  Theorem  5.3  of  [G4]  to  obtain  essentially  the  present 
theorem  without  using  the  artificial,  and  rather  clumsy,  extensions  to  in 
the  proof.  □ 
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Theorem  5.  .2  The  strong  convergence  in  (5.7)  implies  uniform  norm  convergence 
of  the  optimal  feedback  laws: 


HbXpe  -b*hii  0  as  n  — »  ®. 

n  (5.10) 

Proof.  This  follows  from  the  selfadjointness  of  1^  and  PEn  and  the  finite 

dimensionality  of  the  control  space  Rn.  See  equations  (4.23)  and  (4.24)  of 

[Gil.  □ 


Theorem  5 .3  Assume  the  hypotheses  of  Theorem  5.1  but  do  not  assume  (5.6)  or 
(5.9) .  If  ||  IL||  is  bounded  uniformly  in  n,  then  the  Riecati  algebraic 
equation  (3.3)  has  a  nonnegative  selfadjoint  solution  n.  and,  for  each  z  e  E, 
^hPEnz  converges  weakly  to  Hz. 

Proof.  This  is  Theorem  6.7  of  [G3]  ,  whose  proof  is  valid  under  the  hypotheses 
here.  O 

The  main  shortcoming  of  the  weak  convergence  in  Theorem  5.3  is  that  it  does 
not  yield  uniform  norm  convergence  of  the  feedback  control  laws. 


5.2  Convergence  of  the  Approximating  Optimal  Control  Problems  of 


Section  4.2 


For  the  rest  of  this  section,  Aj_,  A,  T(t),  BQ  and  B  are  the  operators 

defined  in  Section  2.  The  operators  An,  Bn,  Qn  and  are  the  operators  in 

the  approximation  scheme  of  Section  4.  In  particular,  II  e  L(E  , E  )  is  the 

n  n 

minimal  nonnegative,  self-adjoint  solution  of  the  Riccati  operator  equation 
(4.26).  According  to  Corollary  4.3  and  Theorem  4.5,  the  Ritz-Galerkin  approx¬ 
imation  scheme  presented  in  Section  4.1  converges  as  required  in  (5.1)  and 
(5.2);  (5.3)  and  (5.4)  follow  from  (4.23)  and  the  definition  Q  =  P„  Q|_  in 
Section  4.2.  Also,  Hypothesis  4.6  guarantees  for  each  n  the  existence  of  the 
required  solution  of  the  Riccati  equation  (5.5)  in  Theorem  5.1. 


Since  11^  is  nonnegative  and  self-adjoint,  its  eigenvalues,  which  are 
also  the  eigenvalues  of  its  matrix  representation,  are  real  and  nonnegative, 
and  its  norm  is  equal  to  its  maximum  eigenvalue. 

Theorem  £.4  If  Q  is  E-coercive  and  dQ  =  0  (i.e. ,  there  is  no  open- loop  damp¬ 
ing),  then  there  is  no  nonnegative  self-adjoint  solution  of  the  Riccati  opera¬ 
tor  equation  (3.3),  and 


I'  HnH  ->  ®  as  n  — »  ®. 


(5.11) 


Proof ♦  Recall  the  operator  A  in  Section  2.1.  By  Theorem  1  of  [G2] ,  there  can 
be  no  compact  operator  C  e  L(E,E)  such  that  X  +  C  generates  a  uniformly 
exponentially  stable  semigroup.  Therefore,  since  a  compact  linear  perturba¬ 
tion  of  X  yields  A,  there  can  be  no  compact  linear  C  such  that  A  +  C  generates 
a  uniformly  exponentially  stable  semigroup. 


Now,  unless  (S.ll)  holds,  there  exists  a  subsequence  such  that  ||  n  || 

J 

is  bounded  in  n . ,  so  that  Theorem  5.3  says  that  there  exists  a  nonnegative 
self-adjoint  solution  II  of  (3.3).  Since  Q  is  coercive.  Theorem  3.3  then 
says  that  the  semigroup  generated  by  A  -  BR  *B*  II  is  uniformly  exponentially 
stable.  But  this  is  impossible  —  BR~*B*  II  is  compact  because  its  rank  is  not 
greater  than  m.  □ 


Theorem  5.5.  Suppose  that  Aq  and  dg(*,*)  are  both  H-coercive.  Then  there 
exist  positive  constants  and  a^,  independent  of  n,  such  that 


-a.  t 

||exp[Antl||  i  ,  t  2  0. 


(5.12) 


Proof.  First,  we  define  A-  and  DA  e  L(V  ,V.)  to  be  the  operators  whose 

un  on  n  n 

matrix  representations  are  M-nKn  and  M-nDn,  respectively.  (See  Section  4.1.) 
The  operator  Ar  is  then 


A  =  0  1  ] 

n  “Aon  "D0nj 


(5.13) 


Since  A^  and  d0  are  H-coercive,  there  exists  a  positive  constant  p, 
independent  of  n,  such  that 


<A0nh'h>H  1  P 1 1 h 1 1 h 


(5.14) 


(5.15) 


for  all  h  e  V  .  Since  dQ  is  continuous  on  V  X  V,  there  exists  a  positive 


constant  y,  independent  of  n,  such  that 


<D0nh*h>H  ^  rllhllj 

n  H  v  (5.16) 

for  al 1  h  e  V  .  The  theorem  follows  then  from  Theorem  2.1.  □ 
n 

Theorem  5. .6.  Suppose  that  Aq  has  an  invariant  subspace  Vq  which  is  a1  to 
invariant  under  the  damping  map  A^Ap  ,  that  eq  =  vq  X  Vq  is  a  stabilizable 
subspace  for  the  control  system,  and  that  the  restrictions  of  AQ  and  d0(»,») 
to  Vq  are  both  H-coercive.  Also,  suppose  that  Vq  has  finite  dimension  nQ  and 
that,  for  each  n  2.  nQ  in  the  approximation  scheme,  the  first  Hq  e^,a  span  VQ 
and  the  rest  are  orthogonal  to  Vq  in  both  V  and  H. 

i)  Then  (3.3)  has  a  nonnegative  solution  II,  and  for  each  n  2.  nQ,  (5.5)  has  a 
nonnegative  self-adjoint  solution  1^.  Also,  Dn  is  bounded  uniformly  in  n,  so 
that  converges  toll  weakly,  as  in  Theorem  5.3. 

ii)  If  Eq  and  E^~  (the  E-orthogonal  complement  of  EQ)  are  invariant  under  Q, 
and  if  the  part  of  the  open-loop  system  on  EQ  is  observable  with  the  measure¬ 
ment  Qz,  then  (5.6)-(5.8)  hold  as  in  Theorem  5.1. 

Proof.  We  will  invoke  Theorem  3.7  to  establish  the  existence  of  the  uniform 
bounds  and  decay  rates  needed  in  Theorem  5.1.  In  the  approximating  optimal 
control  problems,  the  part  of  the  control  system  on  EQ  is  the  same  for  each  n; 
the  approximation  of  the  control  system  takes  place  on  E^.  We  can  write  Er  = 
Eq  e  E^n,  where  E^n  is  the  orthogonal  complement  of  Eq  in  E'n,  and  EQ  and  E^ 
clearly  reduce  the  open- loop  semigroup  for  each  n  2  Hq. 

For  the  part  of  the  open-loop  system  on  E0n>  Theorem  5.5  establishes 


I 
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positive  and  o^»  independent  of  n,  for  (3.17).  Also,  we  have  a  £ 
independent  of  n  for  (3.18)  because  =  PgnB  and  QR  =  PEnQ | En-  Therefore,  i) 
follows  from  Theorem  3.7  i)  and  Theorem  5.3. 

The  definition  of  and  the  requirement  on  where  the  various  basis  vec¬ 
tors  must  lie  imply  that  EQ  and  E^n  reduce  Qn  if  EQ  and  E^  reduce  Q  and  that 
the  restriction  of  Q  to  EQ  is  the  same  for  all  n.  Therefore,  ii>  follows  from 
Theorem  3.7  ii)  and  Theorem  5.1.  □ 

Remark  5. .2.  In  applications,  the  subspace  in  Theorem  5.6  usually  contains 
rigid-body  modes.  The  theorem  includes  the  case  where  both  AQ  and  dQ  are  H- 
coercive  on  all  of  V  (no  rigid-body  modes  and  all  modes  damped).  In  this 
case,  VQ  is  the  trivial  subspace.  O 

Remark  5.8.  Otherwise,  for  applications  to  flexible  structures.  Theorem  5.6 
usually  requires  two  things:  first,  modal  damping  must  be  modeled  for  the 
structure,  so  that  the  natural  modes  remain  uncoupled  in  the  open- loop  system; 
second,  the  natural  mode  shapes  must  be  used  for  the  basis  functions  in  the 
approximating  optimal  control  problems.  Although  these  requirements  may  seem 
restrictive  from  a  mathematical  standpoint,  such  modeling  and  approximation 
predominate  in  engineering  practice.  Also,  we  get  our  strongest  convergence 
results  under  these  conditions.  For  applications  where  the  basis  vectors  are 
not  the  natural  mode  shapes,  the  following  theorem  is  useful.  □ 

Theorem  5.9.  Suppose  that  AQ  +  BqBq  and  dQ  +  BqBq  are  H-coercive.  Then  (3.3) 
has  nonnegative  solution  II,  for  each  n  (5.5)  has  a  nonnegative  self-adjoint 
solution  11^,  and  II  “nil  is  bounded  uniformly  in  n.  Hence  Theorem  5.3 
applies.  Furthermore,  if  Q  is  E-coercive,  then  (5.6)— (5.9)  hold  in  Theorem  5.1. 


Proof.  The  required  bounds  follow  from  Theorem  3.10  and  the  proof  of  Theorem 
|  5.5.  Although  we  took  A ^  =  BqBq  in  Theorem  3.10,  this  is  not  necessary  in  the 

final  result,  since  all  bounded  self-adjoint  operators  on  H  that  make  Aq  + 
A^  coercive  yield  equivalent  norms  for  V.  Q 


Theorem  5.10.  If  (5.7)  holds  for  each  z  e  E,  then 


,,fin  "  f<Mv  0. 


i"V 


(5.17a) 


I  l«in  “  I  lH  “>  0,  as  n  «, 


(5.17b) 


where  f^  and  g^  are  the  functional  gains  in  (3.16),  and  fin  and  gin  are  the 
approximating  functional  gains  in  (4.31)  and  (4.41). 


Proof.  The  result  follows  from  (4.31).  □ 


Note  that  (5.10)  and  (5.17)  are  equivalent 


6.  Example 


6.1  The  Control  System 

One  end  of  the  uniform  Euler- Bernoulli  beam  in  Figure  6.1  is  attached 
rigidly  (cantilevered)  to  a  rigid  hub  (disc)  which  is  free  to  rotate  about  its 
center,  point  0,  which  is  fixed.  Also,  a  point  mass  iS  attached  to  the 
other  end  of  the  beam.  The  control  is  a  torque  u  applied  to  the  disc,  and  all 
motion  is  in  the  plane. 


Figure  6.1.  Control  System 


p  =  hub  radius 
i  =  beam  length 

*0  =  hub  moment  of  inertia  about  axis 
perpendicular  to  page  through  0 

mb  =  beam  mass  per  unit  length 
m^  =  tip  mass 

El  =  product  of  elastic  modulus  and  second  moment 
of  cross  section  for  beam 
fundamental  frequency  of  undamped  structure 


10  in 
100  in 

100  slug  in^ 

.01  slug/in 
1  slug 

13,333  slg  in3/3ec2 
.9672  rad/sec 


Table  6.1  Structural  Data 


The  angle  0  represents  the  rotation  of  the  disc  (the  rigid-body  mode), 
w(t,s)  is  the  elastic  deflection  of  the  beam  from  the  rigid-body  position,  and 

wl<t)  is  the  displacement  of  m^  from  the  rigid-body  position.  For  technical 
reasons,  we  do  not  yet  impose  the  condition  w1(t)  =  w(t,jf);  more  on  this 
later. 


The  control  problem  is  to  stabilize  rigid-body  motions  and  linear  (small) 
transverse  elastic  vibrations  about  the  state  0=0  and  w  =  0.  Our  linear 
model  assumes  not  only  that  the  elastic  deflection  of  the  beam  is  linear  but 
also  that  the  axial  inertial  force  produced  by  the  rigid-body  angular  velocity 
has  negligible  effect  on  the  bending  stiffness  of  the  beam.  The  rigid-body 
angle  need  not  be  small. 

For  this  example,  it  is  a  straight  forward  exercise  to  derive  the  three 
coupled  differential  equations  of  motion  in  0,  w  and  Wj  ,  ancj  they  do  have  the 


« 


« 
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form  (2.1').  However,  to  emphasize  the  fact  that  we  do  not  use  the  explicit 
partial  differential  equations,  we  will  not  write  these  equations  here. 

Rather,  we  will  write  only  what  is  normally  needed  in  applications:  the 
kinetic  and  strain-energy  functionals,  the  damping  functional  and  the  actuator 
influence  operator. 


Remark  2.1  applies  to  this  example,  and  to  most  examples  with  complex 
structures.  The  generalized  displacement  vector  is 


x  =  (©,w,wx)  e  H0  =  R  L2(0,!>X  R. 

The  kinetic  energy  in  th  3  system  is 

•  • 

Kinetic  Energy  =  V2  <x,x>H 
where  H  is  HQ  with  the  inner  product 

>H  “  “b  A) 


A  jf  A 

<x,x>u  -  [w+(r+s)0] [w+(r+s)«lds 


(6.1) 


(6.2) 


(6.3) 


+  +  m1[w1+(r+J?)e][^1+(r+i)^]. 

As  in  most  applications,  we  need  not  write  the  mass  operator  explicitly,  but 
there  exists  a  unique  selfadjoint  linear  operator  Mq  on  HQ  such  that 


<x,1bH  =  <MqX,$>jj  . 


(6.4) 


It  is  easy  to  3ee  that  Mq  is  bounded  and  coercive.  Hence  Hq  and  H  have 
equivalent  norms. 

The  input  operator  for  (2.1')  (which  maps  R  to  HQ)  is 
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(1.0,0). 


(6.5) 


Since  we  multiply  (2.1')  by  to  get  (2.1),  the  input  operator  for  (2.1)  is 


.-1 


(Mq  Bq).  Note  that 

<mo1bo>*H 


V 


(6.6) 


where  (M^Bq)*11  is  the  H-adjoint  of  (Mj^Bg)  and  Bo  is  fche  Ho"adJ°int  of  Bo* 

Remark  2.2  also  applies  here.  The  only  strain  energy  is  in  the  beam  and 
is  given  by 


Strain  Energy  =  V2  a(x,x) 


(6.7) 


with 


a(x,x)  =  El/*  w*  «$' »ds. 


(6.8) 


where  C)  =  3Z( ')/dsz( ') .  To  make  a(*,‘)  into  an  inner  product,  we  must 
account  for  rigid-body  rotation.  Thus  we  set 


<x,x>. 


a(x,x)  + 


0$ 


(6.9) 


and  define 


V  =  (x  =  (0,0, 0(1)):  0  s  H2(0,jf) ,  0(0)  =  0'  (0)  =  0  } 


(6.10) 


Also,  we  have 


<x,$>. 


a(x,x)  +  <B0B*x,$>h 


(6.11) 


=  a(x.x)  +  <(M"1B0)(M^1B0)*Hx,x>H, 
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so  that  Aj  =  BqBq,  or  (M"1^)  (M^1B0> *h,  depending  on  whether  the  HQ  or  the  H- 
inner  product  is  used  in  computing  the  V- inner  product.  But  we  need  neither 
nor  AQ  explicitly.  Mg  need  only  (£.8)  god  (£.9) ,  along  with  (£.3) ,  compute 
ihs  xequlced  inner  product?. 

As  mentioned  in  Remark  2.2,  the  operator  AQ  can  be  defined  now  by  (2.6), 
and  the  stiffness  operator  is  AQ  =  A^  -  Aj.  Using  the  HQ-inner  product  in 
(2.6)  yields  the  AQ  for  (2.1'),  and  using  the  H-inner  product  yields  the  AQ 
for  (2.1),  which  is  M^1  following  the  Aq  for  (2.1').  The  Aq  for  (2.1')  is 
quite  simple,  and  the  reader  might  write  it  out.  We  will  not,  so  that  no  one 
will  think  that  we  use  it.  We  will  point  out  that  D(AQ)  requires  both  the 
geometric  boundary  conditions  in  V  and  the  natural  boundary  condition  w'*(t,JP) 
=  0;  i.e. ,  zero  moment  on  the  right  end. 


Remark  6.1.  That  the  geometric  boundary  conditions 
w(t,0)  =  w'(t.O)  =  0 


and 


(6.12) 


w(t,Jf)  =  wx(t) 


(6.13) 


are  imposed  in  V  but  not  in  H  —  i.e.,  on  the  generalized  displacement  but  not 
on  the  generalized  velocity  —  is  common  in  distributed  models  of  flexible 
structures.  The  natural  norm  for  expressing  the  kinetic  energy  of  distributed 
components  is  the  L^-norm,  which  cannot  preserve  constraints  on  sets  of  zero 
measure.  Because  the  strain  energy  involves  spacial  derivatives,  the  stronger 
strain-energy  norm  can  preserve  the  geometric  boundary  conditions  (although, 
as  for  the  boundary  slope  of  an  elastic  plate,  the  V-norm  may  impose  some  of 
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these  boundary  conditions  in  an  L 2  rather  than  a  pointwise  sense).  The 
strain-energy  norm  is  based  on  the  material  model  of  the  distributed  com¬ 
ponents  of  the  system,  and  it  should  cot  be  surprising  that  such  a  norm  is 
required  to  connect  the  various  structural  components.  Q 


We  assume  that  the  beam  has  Voigt-Kelvin  viscoelastic  damping  [C2] ,  so 
that  the  damping  operator  in  (2.1)  is 


D0  “ 

where  c0  is  a  constant.  This  means  that  the  damping  functional  is 
dQ(x,x)  =  cQ  a(x,x),  x,x  e  V. 


(6.14) 


(6.15) 


6.2  The  Optimal  Control  Problem 


We  take  Q  =  I  in  the  performance  index  in  (3.1).  This  means  that  the 
state  weighting  term  <Qz,z>g  is  twice  the  total  energy  in  the  structure  plus 
the  square  of  the  rigid-body  rotation.  Since  there  is  one  input,  the  oontrol 
weighting  R  is  a  scalar. 


According  to  (3.33),  the  optimal  control  has  the  feedback  form 


u(t )  =  -  <f,x(t)>v  -  <g,x(t)>H 

where  x(t)  has  the  form  (6.1),  and 

f  =  (<xf,0f,pf)  =  R"lniB0  e  V, 


(6.16) 


(6.17a) 


A 


58 


8  =  (WV  =  r-1i^bo8H 


(6.17b) 


Note  that  0^.  =  0j.(Jf)  is  not  used  in  the  control  law — recall  (6.8)  and  (6.9) 


6 .3  Approximation 


Our  approximation  of  the  distributed  model  of  the  structure  is  based  on  a 
finite  element  approximation  of  the  beam  which  uses  Hermite  cubic  splines  as 
basis  functions  ([S1,S4]).  These  are  the  basis  functions  most  commonly  used  in 
engineering  finite  element  approximations  of  beams.  The  splines  and  their 
first  derivatives  are  continuous  at  the  nodes.  Because  the  basis  vectors  ej 
in  the  approximation  scheme  in  Section  4  must  be  in  the  space  V  defined  in 
(6.10),  we  write  them  as 


ex  =  (1,0,0), 


(6.18a) 


®j  =  J  =  2,  3,  .  .  ., 

where  the  0j's  are  the  cubic  splines.  When  we  use  n0 
the  beam,  there  are  2n  linearly  independent  splines. 

C 

body  mode,  the  order  of  approximation  is  n  =  2n  +  1. 

© 


(6.18b) 

elements  to  approximate 
Thus,  with  the  rigid- 


For  the  numerical  solution  to  the  optimal  control  problem,  we  have  only 
to  plug  into  the  formulas  of  Section  4.  The  matrices  in  (4.3)  are  calculated 
according  to  (6.3),  (6.8)  and  (6.9),  with  BQ  given  by  (6.5).  In  particular, 

Kn  =  [a(e,  ,e,)]  ,  Dn  =  c0  Kn,  Mn  =  [<e.,ei>H], 

i  j  0  i  j  H  (6  19) 

Bj  =  [10  0  •**  0]T  =  [<ei,^1(l,0,0)>H]  =  [<e1,(l,0,0)>H  ] 

°(6.20) 
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Note  that  the  first  row  and  column  of  Kn  are  zero.  The  matrix  Kn  in  (4.32)  is 
Kn  with  1  added  to  the  first  element.  The  matrices  An  and  Bn  are  given  by 
(4.5)  and,  sinoe  Q  =  I,  the  matrix  Qn  is  the  Wn  in  (4.33).  With  these 
matrices,  we  solve  the  Riccati  equation  (4.40)  and  use  (4.41)  and  (4.45)  to 
compute  the  approximations  to  the  functional  gains,  which  are 


fn  =  (afn’^fn'pfn) * 


(6.21a) 


(o_  #0  iP__) . 

n  gn  gn  gn 


(6.21b) 


For  convergence,  we  satisfy  all  the  hypotheses  of  Theorem  5.9.  In  par¬ 
ticular,  since  Q  is  the  identity  on  E,  it  is  coercive.  Theorem  5.9  implies 
that  the  solutions  to  the  finite  dimensional  Riccati  equations  converge  as  in 
Theorem  5.1  and  that  the  functional  control  gains  converge  as  in  theorem  5.10. 


Remark  £.2.  It  might  appear  that  the  hypotheses  of  Theorem  5.6  hold  with  = 

1,  but  not  so.  For  j  1  2,  is  orthogonal  to  ex  in  HQ  and  V  not  in  H. 

Recall  (6.1)-(6.3) ,  (6.9)-(6.11)  and  (6.18) .  O 


6.4  Numerical  Results 

Figures  6.2a  and  6.2b  shew  the  computed  functional  gain  kernals  4^’  and 
’  6^  for  the  damping  coefficient  c0  =  10-4*  the  control  weighting  R  =  1,  and 

ne  -  2,  3,  4,  5  and  8  beam  elements.  Table  6.2  lists  the  corresponding  scalar 

components  of  the  gains.  For  oQ  -  10-4  and  R  =  05,  the  convergence  Is  slower, 

f  as  discussed  below.  To  show  the  complete  story  of  convergence.  Figures  6.3a 

and  6.3b  and  Table  6.3  show  the  results  for  n  =  2, 3  ,4, 5, 8,  and  Figures  6.4a 

and  6.4b  and  Table  6.4  shew  the  results  for  n  =  4,6,8,10. 

G 

9  9 

We  have  plotted  0fn  because  the  second  derivative  appears  in  the 
strain-energy  inner  product  in  (6.8)  and  (6.9)  and  converges  in  H^O.Jf). 

Note  that,  since  the  Hermlte  cubic  splines  have  discontinuous  second  deriva- 

9  9 

tlves  at  the  nodes,  the  approximations  to  0f  are  discontinuous  at  the  nodes. 

Although  H2-convergence  guarantees  only  ^-convergence  for  0^.  it  can  be 

shown  that  <k' '  converges  uniformly  on  [0,|]  for  this  problem, 
in 

The  tables  omit  pfn  to  emphasize  the  fact  that  it  does  not  appear  in  the 
feedback  law  and  the  fact  that  the  convergence  of  pfn  is  not  an  independent 
piece  of  Information  about  the  convergence  of  the  control  gains;  since  0fn(°) 

=  0fn(O)  =  0,  the  convergence  of  0fn  implies  the  convergence  of  PfR  =  0fn  (j() . 
On  the  other  hand,  although  PgR  =  0^  (J?)  for  each  n,  the  H-norm  convergence 
of  gR  does  not  enforce  this  condition  in  the  limit,  as  the  V-norm  convergence 
of  fR  enforces  .  Hence,  as  far  as  we  can  tell  from  our  results  in 

Sections  3.5,  p^n  is  an  independent  indicator  of  the  convergence  of  the  con¬ 
trol  gains,  as  well  as  being  used  in  the  control  law  in  (6.16).  However,  the 
behavior  of  0^  in  Figures  6.2b,  6.3b  and  6.4b  suggests  that  gn  converges  in 
V.  Stronger  results  on  the  continuity  of  0g  and  the  convergence  of  0gn(j() 
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should  follow  from  a  theorem  stating  that,  because  the  open- loop  semigroup 
generator  A  is  analytic,  the  solution  to  the  infinite  dimensional  Riccati 
equation  maps  all  of  E  into  D(A*) .  The  fact  that  $^(J?)  converges  to  zero  in 
Figure  6.2a  also  suggests  such  a  theorem,  but  we  have  not  proved  it. 

With  the  state  weighting  Q  fixed,  the  two  factors  that  determine  the  rate 
of  convergence  are  c^  and  R.  Although  we  have  used  splines  to  approximate  the 
beam,  the  relation  between  the  convergence  rate  and  Cq  and  R  probably  can  be 
interpreted  best  in  terms  of  the  number  of  natural  modes  of  the  structure  that 
the  optimal  infinite  dimensional  controller  really  controls.  Strictly  speak¬ 
ing,  the  controller  controls  all  modes,  but  the  functional  gains  lie  essen¬ 
tially  in  the  span  of  some  finite  number  of  modes.  This  would  be  the  number 
of  modes  required  for  convergence  of  the  gains  if  we  used  the  natural  modes  as 
the  basis  vectors  in  the  approximation.  The  rest  of  the  modes  are  practically 
(but  not  exactly)  orthogonal  to  the  functional  gains,  so  that  the  optimal 
feedback  law  essentially  ignores  them.  In  general,  the  lighter  the  damping, 
the  more  modes  that  will  be  controlled  for  given  Q  and  R;  the  cheaper  the  con¬ 
trol,  the  more  modes  that  will  be  controlled  for  given  Q  and  cQ.  The  question 
of  the  convergence  of  the  finite  element  approximation  to  the  functional  gains 
becomes  then  a  question  of  how  many  modes  the  optimal  control  law  really  wants 
and  how  many  elements  it  takes  to  approximate  those  modes. 

Numerical  experience  with  optimal  control  of  flexible  structures  has 
shown  this  modal  interpretation  of  the  convergence  of  the  approximating  con¬ 
trol  laws  to  be  very  useful,  and  that  it  is  difficult  to  improve  upon  the 
natural  modes  as  basis  vectors  for  the  approximation  scheme  (see  [G5lX  How¬ 
ever,  whether  the  natural  modes  are  always  or  almost  always  the  best  basis 


( 
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vectors  is  an  open  question.  We  use  the  cubic  splines  here  to  demonstrate 
that  a  standard  finite  element  approximation  works  quite  well.  Also,  to  use 
the  natural  modes  as  basis  vectors  here,  we  first  would  have  to  compute  them 
using  a  finite  element  approximation  —  as  in  most  real  problems  —  and  we  do 
not  know  in  advance  which  or  hew  many  modes  are  needed.  On  the  other  hand,  if 
the  most  important  natural  modes  are  determined  from  experiment,  then  modal 
approximation  should  be  best. 
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ne 

°fn 

a__ 

gn 

Pgn 

2 

1.000 

.2141 

-22.459 

3 

1.000 

.2396 

-25.221 

4 

1.000 

.2496 

-26.331 

5 

1.000 

.2534 

-26.786 

8 

1.000 

.2561 

-27  .041 

Xafels  6.2 

.  Scalar  Components 

of  Functional  Control  Gains 

Damping  coefficient 

°o  - 

;  control  weighting  R  =  1 

number  of 

elements 

ng  =  2 ,  3 , 

4.  5,  8 

ne 

“fn 

agn 

^gn 

2 

4.4721 

1.0136 

-108.27 

3 

4.4721 

1.1770 

-126  .40 

4 

4.4721 

1.2440 

-133 .87 

5 

4.4721 

1.2781 

-137.57 

8 

4.4721 

1.3106 

-141.15 

Table  6 .3 , 

.  Scalar  Components 

of  Functional  Control  Gains 

Damping  coefficient 

cQ  =  io~4; 

oontrol  weighting  R  =  .05 

number  of 

elements 

ne  =  2.  3, 

4.  5,  8 

n 

e 

°fn 

ag  n 

^gn 

4 

4.4721 

1.2440 

-133.87 

6 

4.4721 

1.2973 

-139.69 

8 

4.4721 

1.3106 

-141.15 

10 

4.4721 

1.3141 

-141.54 

Table  6 .4 . 

,  Scalar  Components 

of  Functional  Control  Gains 

Damping  coefficient 

co  = 

control  weighting  R  =  .05 

number  of 

elements 

ne  =  4.  6, 

8,  10 
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Figures  6.5a  and  6.5b  and  Table  6.5  represent  attempts  to  compute  an 
optimal  control  law  for  the  structure  when  R  =  .05  but  eQ  =  0.  Since  Q  is  the 
identity  operator  in  E  and  hence  ooercive,  Theorem  5.4  says  that  no  optimal 
control  law  exists  and  that  the  norm  of  the  solution  to  the  finite  dimensional 
Riccati  equation  grows  without  bound  as  the  number  of  elements  increases. 

This  is  reflected  in  the  nonconvergence  of  agn,  0^  and  0gn>  although  afn  con- 

9  $ 

verges  and  the  convergence  of  is  unclear. 

In  applications  where  the  structural  damping  is  not  known,  except  that  it 
is  very  light,  it  is  tempting  and  not  uncommon  engineering  practice  to  assume 
zero  damping  in  the  design  of  a  control  law  for  the  first  few  modes,  while 
trusting  whatever  damping  is  in  the  higher  modes  to  take  care  of  them.  How¬ 
ever,  if  high  performance  requirements  (large  Q)  or  ooupling  between  modes  in 
the  closed-loop  system  necessitate  a  control  law  based  on  a  more  accurate 
approximation  of  the  structure.  Theorem  5.4  and  the  current  example  warn  that 
the  higher-order  control  laws  are  likely  meaningless  and  rather  strange  if  no 
damping  is  modeled. 

We  should  note  that  we  have  seen  similar  problems  [G9]  where  remains 
bounded  and  the  gains  converge  for  zero  damping  but  finite-rank  Q.  In  such 
cases.  Theorem  5.3  says  that  an  optimal  control  law  exists  for  the  distributed 
model  of  the  structure  and  that  the  finite  dimensional  control  laws  converge 
to  an  optimal  infinite  dimensional  control  law.  Also,  Balakrishnan  [B2]  has 
shown  that  an  infinite  dimensional  optimal  control  law  exists  for  no  damping 
when  Q  =  BB*. 
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ne 

°fn 

°gn 

^gn 

2 

4.4721 

1.0516 

-112.23 

3 

4.4721 

1.3061 

-140.18 

4 

4.4721 

1.4758 

-159.11 

5 

4.4721 

1.5996 

-172.64 

8 

4.4721 

1.8407 

-199.39 

Table  £.£.  Scalar  Components  of  Functional  Control  Gains 
Zero  damping;  control  weighting  R  =  .05 
number  of  elements  ng  =  2,  3,  4,  5,  8 
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7 .  Iks  optimal  Infinite  Dimensional  Estimator.  C<?BP<jyi9ftfcor  and  Closed-loop 
System 

As  in  Sections  3  and  5,  we  will  state  some  initial  definitions  and 
results  for  an  arbitrary  linear  control  system  on  a  Hilbert  space  in  Subsec¬ 
tion  7.1,  and  then  discuss  implications  for  flexible- structure  control  in  Sub¬ 
section  7 .2 . 


7.1  The  Generic  Problem 


Let  A,  T(t)  and  B  be  as  in  Subsection  3.1,  with  E  an  arbitrary  real  Hil¬ 
bert  space.  The  differential  equation  corresponding  to  (3.2)  is,  of  course. 


z(t)  =  Az(t)  +  Bu(t),  t  >  0. 


(7.1) 


We  assume  that  we  have  a  p-dimensional  measurement  vector  y(t)  given  by 


y(t)  =  C  u(t)  +  Cz ( t ) , 

0  (7.2) 

where  CQ  e  L(Rm.  Rp)  and  C  e  L(E,RP)  for  some  positive  integer  p. 

Definition  1.1.  For  any  F  e  L(RP,E),  the  system 

z(t)  =  Az(t)  +  Bu(t)  +  F[y(t)  -  C0u(t)-Cz(t)l ,  t  >  0, 

V  (7.3) 


will  be  called  an  observer,  estimator  (we  use  the  terms  interchangeably) ,  for 

A  A 

the  system  (7.1)— (7.2) .  Let  S(t)  be  the  semigroup  generated  by  A-FC.  The 

A 

observer  in  (7.3)  is  strongly  (uniformly  exponentially)  stable  if  S(t)  is 
strongly  (uniformly  exponentially)  stable.  □ 


To  justify  this  definition,  we  write 
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e(t) 


z(t)  -  2(t) 


(7.4) 


and,  with  (7.1)-(7.3),  obtain 

e(t)  =  §(t)e(0),  t  1  0. 


(7.5) 


Of  course,  an  observer,  or  estimator,  is  necessary  because  the  full  state 
z(t)  will  not  be  available  for  direct  feedback,  and  the  feedback  control  must 
be  based  on  an  estimate  of  z(t).  When,  as  in  this  paper,  the  desired  control 
law  has  the  form 


u(t)  =  -Fz(t) 


(7.6) 


for  some  F  e  L(E,lP),  the  observer  in  (7.3)  can  be  used  to  construct  z(t)  from 
the  measurement  in  (7.2)  and  then  the  control  law  in  (7.6)  can  be  applied  to 
z(t).  The  control  applied  to  the  system  is  then 


u(t)  =  -Fz ( t ) , 


(7.7) 


and  the  resulting  closed-loop  system  is 


t  l  0, 


(7.8) 


where  S^.^tt)  is  the  semigroup  generated  on  E  x  E  by  the  operator 


A  -BF  ] 

[A-BF-$C]|*  D*Aa»'<»)  =  D(A)XD(A). 


(7.9) 


With  the  estimator  error  e(t)  defined  by  (7.4),  it  is  easy  to  show  that 

(7.8)  is  Equivalent  to  (7.5)  and 
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Z(t)  =  (A-BP)z(t)  +  BFe (t ) ,  t  >  0, 


(7.10) 


where  (A-BF)  generates  a  semigroup  S(t)  on  E.  Also,  it  is  easy  to  prove  the 
following. 

Theorem  7 .2  ♦  Suppose  that  there  exist  positive  constants  Mj,  cij  and  02 
such  that 

-a.  t 

IIS(t)||  i  ^  e 


I  |S(t) 1 1  i  Mje 


t  1  0. 


(7.11) 


Then,  for  each  real  <  mintaj,  aj) ,  there  exists  a  constant  such  that 


II  S„,„(t)|  |  1  Mje  ,  t  2  0. 


(7.12) 


Also, 


<*(A«,J  =  o(A-BF)  U  o(A-FC), 

where  ©(A^,^)  is  the  spectrim  of  A^,^.  □ 


(7.13) 


The  observer  in  (7.3)  and  the  control  law  in  (7.7)  constitute  a  compensa¬ 
tor  for  the  oontrol  system  in  (7.1)  and  (7.2).  The  transfer  function  of  this 
compensator  is 


«(s) 


=  -F  (si- [A-BF 


A  _1  A 

+  F(C  F-C)J)  XF 


(7.14) 


which  is  an  m  x  p  matrix  function  of  the  complex  variable  s.  When  E  has 
infinite  dimension,  the  compensator  transfer  function  is  irrational,  except  in 
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degenerate,  usually  unimportant  cases 


The  foregoing  definitions  of  this  section  and  Theorem  7.1  are  straight¬ 
forward  generalizations  to  infinite  dimensions  of  observer-controller  results 
in  finite  dimensions.  Balas  [B3]  and  Schumacher  [S2]  have  used  similar  exten¬ 
sions. 


Now  suppose  that  F  is  chosen  as 


t  -  ftc’S-1, 

(7.15) 

A 

where  II  e  L(E,E)  is  the  minimal  nonnegative  selfadjoint  solution  to  the  Ric- 
catl  equation 


A 

AH 


A 

n 


A 

-  n 


A- 1  A 

r  cn 


A 

Q 


=  0, 


(7.16) 


with  Q  e  L(E, E)  nonnegative  and  selfadjoint  and  R  e  L(RP,RP)  symmetric  and 
positive  definite.  Theorem  3.3  (with  A,  B,  Q,  R,  II  and  S(t)  replaced  by 
A*. C*,  Q,  R,  6.  and  S*(t))  gives  sufficient  conditions  f or  ft  to  exist  and  for 

A*  A 

the  semigroup  S  (t)  —  and  equivalently  its  adjoint,  the  S(t)  generated  by 

A  *  A— J 

A-  D.  C  R  c  —  to  be  uniformly  exponentially  stable. 


Definition  2.3.  When  the  control  gain  operator  is 
F  =  R-1B*II  . 


(7.17) 


with  II  the  solution  to  the  Riccati  equation  (3.3),  and  the  observer  gain 
operator  is  given  by  (7.15)  and  (7.16),  we  will  call  the  compensator  consist¬ 
ing  of  the  observer  in  (7.3)  and  the  control  law  in  (7.7)  the  optimal  infinite 
dimensional  compensator,  and  (7.8)  the  optimal  closed-loop  system.  □ 


Ik 


*  Az  +  Bu 


CqU  +  Cz 


Control  System 


[A-BF-F (CQ  F-C)]  z  +  F  y 


Optimal  Infinite  Dimensional  Compensator 


Figure  7.1.  Optimal  Closed-loop  System 


Remark  7 .4 .  The  infinite  dimensional  observer  defined  by  (7.3),  (7.15)  and 
(7.16)  is  the  optimal  estimator  for  the  stochastic  version  of  (7.1)  and  (7.2) 

when  (7.1)  is  disturbed  by  a  stationary  gaussian  white  noise  process  with  zero 

A 

mean  and  covariance  operator  Q  and  the  measurement  in  (7.2)  is  contaminated  by 

A 

similar  noise  with  covariance  R.  For  infinite  dimensional  stochastic  estima¬ 
tion  and  control,  see  [Bl,  C4] .  When  the  state  weighting  operator  Q  in  (3.1) 
is  trace  class,  the  optimal  infinite  dimensional  compensator  minimizes  the 
time-average  of  the  expected  steady-state  value  of  the  integrand  in  (3.1). 
Existing  theory  for  stochastic  control  of  infinite  dimensional  systems 
requires  trace-class  Q,  but  we  have  a  well  defined  compensator  for  any  bounded 
nonnegative  selfadjoint  Q  and  ft,  as  long  as  the  solutions  to  the  Riccati 


equations  exist.  As  the  next  two  sections  show  (without  assuming  trace-class 

A 

Q) ,  the  infinite  dimensional  compensator  is  the  limit  of  a  sequence  of  finite 
dimensional  compensators,  each  of  which  can  be  interpreted  as  an  optimal  LQG 
compensator  for  a  finite  dimensional  model  of  the  structure.  Therefore,  we  do 
not  require  trace-class  Q  in  our  definition  of  the  optimal  compensator,  even 
though  this  compensator  solves  a  precise  optimization  problem  only  when  Q  is 
trace  class. 

This  paper  is  concerned  primarily  with  hew  the  finite  dimensional  compen¬ 
sators  oonverge  to  the  infinite  dimensional  compensator,  and  the  analysis  of 
this  convergence  requires  only  the  theory  of  infinite  dimensional  Riccati 
equations  for  deterministic  optimal  oontrol  problems  and  the  corresponding 
approximation  theory.  While  the  stochastic  interpretation  of  the  infinite 
dimensional  compensator  and,  in  Section  8.2,  of  the  finite  dimensional  estima¬ 
tors  should  be  motivational,  nothing  in  the  rest  of  the  paper  depends  on  a 

A  A 

stochastic  formulation.  We  assume  that  the  operators  Q,  R,  Q  and  R  are  deter¬ 
mined  by  3ome  design  criteria.  In  many  engineering  applications,  determinis¬ 
tic  criteria  such  as  the  stability  margin  and  robustness  of  the  closed-loop 
system,  rather  than  a  stochastic  performance  index  and  an  assumed  noise  model, 

A  A  p| 

govern  the  choice  of  Q,  R,  Q  and  R.LI 
7.2  Application  to  Structures 

For  the  rest  of  the  paper,  E  =  V  X  H  as  in  Section  2,  and  A  and  B  are  the 
operators  defined  there. 

The  measurement  operator  C  in  (7.2)  now  must  have  the  form 


tci  CJ 


(7.18) 


where  e  L(V,RP)  and  e  L(H,RP).  Hence,  if  we  denote  by  (C(x,x))^  the  i 


th 


component  of  the  p-vector  C(x,x),  for  (x,x)  e  E,  then  there  must  exist  cJi  e  V 
and  e2i  e  H  such  that 


(C(x#x))^  *  ^C^»x)y  +  ^2i'  X^H'  i  ~  1, 


(7.19) 


Also,  the  estimator  gain  operator  F  is  given  by 


A  P  A  A 

-  I  (fVgJy, 

i=l  1  1  1 


(7.20) 

A 


for  y  -  [yj  y2...  y  ]  e  Rp,  where  the  functional  estimator  gains  f±  and  g^ 
are  elements  of  V  and  H,  respectively. 


For  the  optimal  estimator  gains,  we  can  partition  II  as 


fi  - 


-A  A  _ 

*0  V 

A*  A 

«i  ^ 


and  use  (7.15)  and  (7.19)  to  get 


(7.21) 


A 

fi  " 


P  A_ I  A  A 

j?i  R  +  D^le2j)  * 


(7.22a) 


A 


P  A_i  A*  A 

1  =  j?i  (R  )ij(IIlClJ  +  I^2c2j)  ’  1  =  l’2 


(7.22b) 


Now  let  us  partition  Q  as  in  (4.34): 
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(7.23) 


a 


In  the  optimal  control  problem,  we  almost  always  have  a  nonzero  QQ  because 
this  operator  penalizes  the  generalized  displacement.  For  the  results  in  this 

A 

paper,  QQ  can  be  nonzero  in  the  observer  problem,  and,  as  in  the  control  prob¬ 
lem,  some  of  the  strongest  convergence  results  for  finite  dimensional  approxi¬ 
mations  can  be  proved  only  for  coercive  However,  if  the  observer  is  to  be 

A 

thought  of  as  an  optimal  filter,  then  Q  should  be  the  covariance  operator  of 
the  noise  that  disturbs  (2.1).  In  this  case,  =  o  and  ^  =  0. 


i 

I 


I 
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8.  Approximation  of  the  Infinite  Dimensional  Estimator 


f 

8.1  The  Approximating  Finite  Dimensional  Estimators 

Here,  the  scheme  for  the  approximation  of  the  flexible  structure  is  that 
in  Section  4.  We  will  construct  on  the  subspace  En  an  estimator  that  approxi¬ 
mates  the  optimal  infinite  dimensional  estimator  of  Section  7,  and  this  esti- 
mator  will  produce  an  n-order  estimate  =  ($n»$n)  of  the  infinite  dimen- 
sional  state  vector  z  =  (x,x).  In  Section  9,  the  the  n-order  compensator 

t*  h 

that  results  from  applying  the  n  n  approximation  to  the  optimal  control  law 
(in  Section  4.2)  to  z  will  approximate  the  optimal  infinite  dimensional  com¬ 
pensator  of  Section  7. 


Hypothesis  8 .1 .  There  exist  a  sequence  Cn  e  L(En»RP)  such  that 
llcnpEn~cll  0  as  n  ->  ® 

A  A*  A 

and  a  sequence  Q  e  L(E  )  ,  Q  =  Q  >  0,  such  that 

n  n  n  n  — 

A  A 

QnPEn  ®  strongly  as  n  — >  00  .□ 


(8.1) 


(8.2) 


*  * 

Hypothesis  8 .2 .  For  each  n,  the  system  (A^.C^)  is  stabil  izable.  In  particu¬ 
lar,  any  unstable  modes  of  the  system  <cn*An)  are  observable.  □ 


The  nth  observer,  or  nth  estimator.  Is 


A  z  B  u  +  F  ( y-Cnu-C  z  ) 

nn  n  n  unn 


where  the  estimator  gain  is 


(8.3) 


A  /V  *A_1 

Fn  -  W  1  . 


(8.4) 


and  is  the  nonnegative  selfadjoint  solution  to  the  Riceati  operator  equa- 


AA  +  *  3n  ■  0-  (8  5) 

Hypothesis  8.2  implies  that  suoh  a  solution  exists  and  is  unique. 


This  representation  of  the  ntn  estimator  as  a  system  on  E^,  with  the 
estimator  gain  determined  by  the  solution  to  a  Riccati  operator  equation,  is 
necessary  for  showing  how  the  sequence  of  finite  dimensional  estimators 
approximate  the  infinite  dimensional  estimator.  However,  on-line  computations 
will  be  based  on  the  equivalent  differential  equation 


n  =  An  ^  +  Bnu  +  F^y-C^u-C1^) 


,2n  .n 


(8.6) 


where  q(t)  c  R  .A  and  B  are  the  matrix  representations  of  the  operators  Af 
and  Bn,  as  in  Section  4,  and  Cn  is  the  matrix  representation  of  Cn. 

A 

The  2nxp  gain  matrix  F  is 


Fn  =  n11  w~n  (C")^-1. 


where  Wn  is  the  2n  X  2n  grammian  matrix  in  (4.33)  and  11°  satisfies 

An n°  +  (An)T  Hn  +  n0  wn(cn)TR_1cnnn  +  on  =  o, 

with  Qn  the  matrix  representation  of  The  relationship  between  zn  = 

(x  ,1c  )  md  is,  of  course, 
n  n 


(8.7) 


(8.8) 
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(8.9) 


x  (t)  = 

n' 


n  a 

I  5 

i=l 


and 

•?  -  t?T  |TiT. 

(8.10) 

Since  the  matrix  representations  of  Ar  and  A*  are  An  and  W-n(An)T  Wn,  respec¬ 
tively,  and  the  matrix  representation  of  C*  is  W_n(Cn)T,  (8.7)  is  the  matrix 

representation  of  (8.4),  and  the  2nX2n  Riccati  matrix  equation  (8.8)  is  the 

A-  A 

matrix  representation  of  (8.5),  with  IP1  the  matrix  representation  of  1^. 
(Recall  that  W~n  is  the  inverse  of  Wn. ) 

As  in  the  control  problem,  we  do  not  solve  the  matrix  representation  of 
the  nth  Riccati  operator  equation  directly  because  the  matrix  representation 
of  a  selfadjoint  operator  in  general  is  not  symmetric.  In  the  duality  between 
the  optimal  control  and  estimator  problems,  (8.5)  and  (8.8)  correspond  to 
(4.26)  and  (4.38),  respectively.  In  (4.39),  we  defined  the  symmetric  matrix 
IT1  =  Wn  tf1  and  then  obtained  the  Riccati  equation  (4.40)  to  solve  for  IP1  . 

We  proceed  in  a  similar  fashion  here,  but  with  an  interesting  difference. 

A  A  /Vj 

Since  IJ^  and  Qr  are  nonnegative  selfadjoint  operators  on  and  IT  and 
Qn  are  their  matrix  representations,  the  matrices  Wn  fi11  and  are  nonnega¬ 

tive  and  symmetric.  Hence,  the  matrices 


=  ihrn 


and 


(8.11) 
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Qn  =  frVn 

(8.12) 

are  nonnegative  and  symmetric.  Substituting  (8.11)  and  (8.12)  into  (8.7)  and 
(8.8)  yields 

Pn  =  ^(C0)7  R-1 

(8.13) 


and 


An§n  +  ^(An}T  _  +  %n  =  0> 

(8.14) 

the  Riccati  matrix  equation  to  be  solved  numerically  in  the  nth  approximation 
to  the  infinite  dimensional  estimator.  In  view  oi  the  relationship  between 
(8.5)  and  (8.8)  and  the  relationship  between  (8.8)  and  (8.14),  we  see  that 
Hypothesis  8.1  guarantees  the  existence  of  a  unique  nonnegative  symmetric 
solution  to  (8.14). 


To  see  the  relationship  between  the  matrices  in  (8.14)  and  the  operators 

In  (8.5)  more  clearly  —  and  the  difference  between  the  current  approximation 

scheme  and  that  used  in  Section  4.2  for  the  oontrol  problem  —  suppose  that  we 
a  A  5!. 

take  Qn  =  PEnQ| En.  Let  Q  be  defined  as  in  (4.36)  and  (4.37)  with  QQ ,  ,  and 

A  A  A 

Q2  replaced  by  Q Q- ,  and  Q2 .  Then 


W"n 


Q 


n 


W  n. 


(8.15) 


For  example,  if  Q  in  the  control  problem  and  3  in  the  estimator  problem  are 
both  equal  to  the  identity,  then  the  Qn  in  (4.35)  -  (4.42)  is  Wn  anJ 

A/ 

Qn  =  W  n.  This  may  seem  suspicious,  but  Subsection  8.2  should  demonstrate 
that  we  are  solving  the  appropriate  estimator  problem  here. 


The  only  thing  missing  now  for  numerical  implementation  of  the  n^ 

1 


82 


r 


estimator,  or  observer,  is  to  give  Cn,  the  matrix  representation  of  C  ,  expli¬ 
citly.  We  write 


cn  =  [cj  c£] 


(8.16) 


where  the  px  n  matrices  and  c£  are,  respectively,  the  matrix  representa¬ 
tions  of  the  operators  and  C2  in  (7.18).  We  can  cover  virtually  all  appli¬ 
cations  by  assuming  C  =  Cl-  ,  in  which  case  the  i^1  column  of  cf  is  the  p- 

n  nn  J- 

y. 

vector  equal  to  C-je^,  and  the  i  column  of  is  the  p-vector  equal  to  C^e^. 

w§  H2M  have  Hie  .afifflBl&fcfi  .set  of  equations  far  jmneciaal  laBleafintafeigB  of 
the  ntlj  state  estimator:  For  online  computation,  the  nth  estimator,  or 

A 

observer,  is  (8.6);  the  gain  matrix  Fn  is  given  by  (8.13)  and  the  solution  to 
the  Riccati  matrix  equation  (8.14).  The  matrices  Qn  and  Cn  are  defined  as 
above. 


8.2  Stochastic  Interpretation  of  the  Approximating  Estimators 

As  we  have  said,  our  approximation  theory  for  the  optimal  estimator  is 
based  on  approximation  of  the  infinite  dimensional  Riccati  equation,  whose 
structure  is  the  same  for  both  control  and  estimator  problems,  and  the  sto¬ 
chastic  properties  of  the  optimal  estimator  problem  never  enter  our  approxima¬ 
tion  theory.  Furthermore,  using  only  the  deterministic  setting  above,  we  will 
proceed,  subsequently,  to  analyze  the  finite  dimensional  estimators  and  the 
compensators  based  upon  them.  Nontheless,  we  should  consider  momentarily  the 
sequence  of  finite  dimensional  stochastic  estimation  problems  whose  solutions 
are  given  by  the  equations  of  the  preceding  subsection. 

First,  recall  how  the  covariance  operator  of  a  Hilbert  space-valued  ran- 
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dom  variable  is  defined.  The  covariance  operator  of  an  E-valued  random  vari¬ 
able  u  is  the  operator  Q  for  which 


expected  value  {<z,<o  >_<£,<«>>-}  =  <Qz/z>.,,  z,  z  e  E. 

E  E  E  (8.17) 

(See  [Bl,  C4] .) 

A 

With  Fn  given  by  (8.4)  and  (8.5),  (8.3)  is  the  Kalman-Bucy  filter  for  the 
system 


z 

n 


A  z 
n  n 


V 


(0 


'n* 


(8.18) 


y  =  C0U  +  Cnzn  +  “0* 


(8.19) 


where  a>n(t)  is  an  En~valued  white  noise  process  with  covariance  operator  Qn 
and  o>0(t)  is  an  Revalued  white  noise  process  with  covariance  operator 
(matrix)  R.  Next,  careful'  inspection  will  show  that  the  filter  defined  by 
(8.6),  (8.13)  and  (8.14)  is  the  matrix  representation  of  the  filter  defined  by 
(8.3),  (8.4)  and  (8.5). 


With  zr  and  q  "elated  as  in  (4.1)  and  (4.4),  (8.18)  and  (8.19)  are 
equivalent  to  the  system 


t\  =  +  Bnu  +  (/, 


(8.20) 


y  *  Cftu  +  Cl  +  <on , 


-,2n 


(8.21) 


where  (/( t)  is  the  R  -valued  noise  process  related  to  “n(t)  by 
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(8.22) 


n 

«>n(t)  =  I  (O^Ct )e±,  </i+n<t)ei). 

Certainly,  a  Kalman-Bucy  filter  for  (8.20)  and  (8.21)  has  the  form  (8.6)  with 
the  filter  gain  given  by  (8.13)  and  (8.14).  This  particular  filter  is  the 
matrix  representation  of  the  filter  defined  by  (8.3),  (8.4)  and  (8.5)  if  and 
only  if  the  matrix  Q1  defined  by  (8.12)  is  the  covariance  of  the  process  (Jit). 

An  A 

Since  Q  is  the  matrix  representation  of  Qn,  straightforward  calculation  using 
(8.12)  and  (8.17)  shows  that  the  Qn  in  (8.12)  is  indeed  the  correct  covariance 
matrix. 

Of  course,  if  wn(t)  and  (/( t )  represent  a  physical  disturbance  to  the 

(2) 

structure,  then  wn(t)  must  have  the  form  (0,&/  (t))  and  the  first  n  elements 

of  (/( t )  must  be  zero,  but  this  is  not  necessary  for  our  analysis. 

The  finite  dimensional  observers  can  be  interpreted  now  as  a  sequence  of 
filters  designed  for  the  sequence  of  finite  dimensional  approximations  to  the 
flexible  structure,  with  the  n1"*1  approximate  system  disturbed  by  the  noise 

A 

process  o>n(t),  whose  covariance  operator  is  Qn<  By  Hypothesis  8.1,  these 

A 

covariance  operators  converge  to  the  operator  Q  of  Section  7.  If  we  have  a 
reliable  model  of  a  stationary,  zero-mean  gaussian  disturbance  for  the  struc¬ 
ture,  then  we  can  take  the  covariance  operator  for  this  disturbance  to  be  $ 
and  think  of  the  infinite  dimensional  observer  as  the  optimal  estimator.  But, 
again,  this  interpretation  is  not  necessary  for  the  rest  of  our  analysis. 

8.3  The  Approximating  Functional  Estimator  Gains 

The  n  estimator  gain  operator  in  (8.4)  has  the  same  form  as  the  infin¬ 
ite  dimensional  estimator  gain  in  (7.15)  and  (7.20).  We  have 
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(8.23) 


for  y  =  [y,  y,  ...  y  e  Rp,  where  the  functional  estimator  gains  £ .  and 
1  *  p  m 

®in  are  elements  =  Hn*  The  matrix  Fn  in  (8.7)  and  (8.13)  is  the  matrix 


representation  of  Fn,  which  means  that,  if  we  write 


F  = 


f,  f„ 

^  1  z  . .  •  p 

p  p  p 

Lp  p  p p. 


where  the  columns  p  ,  p  E  R  ,  then 


(8.24) 


1  PJ  6j* 
j=l  J  J 


—  l,.«.,p. 


(8.25a) 


n  g, 

I  P,  e  i  =  l,...,p. 

j=l  J 


(8.25b) 


For  convergence  analysis,  it  is  useful  to  note  that  fin  and  gin  are  given 
also  by  equations  corresponding  to  (7.22).  With  the  measurement  operator  C 
written  as  in  (7.19)  and  =  c|£  ,  we  have 


JE^R  )ij(  nonPVnCij  +  IIlnPHnC2j)  ' 


(8.26a) 


where 


JJ1(R  )  ij (  ninPVnci j  +  n2nPHnc2j)' 


(8.26b) 


A 

^  = 


A  A  _ 

^On  Dln 

°l’„  in 


(8.27) 
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8.4  Convergence 


Now  we  will  indicate  the  sense  in  which  the  finite  dimensional 
estimators/observers  approximate  the  infinite  dimensional  estimator  in  Section 
7.  As  we  have  said,  implementation  of  the  n^h  estimator  is  based  on  (8.6), 
(8.13)  and  (8.14),  but  convergence  analysis  is  based  on  the  equivalent  system 
(8.3),  (8.4)  and  (8.5).  The  question  then  is  how  the  observer  in  (8.3),  with 
gain  given  by  (8.4)  and  (8.5),  converges  to  the  observer  in  (7.3)  with  gain 
given  by  (7.15)  and  (7.16). 


Recall  Hypothesis  8.1,  and  recall  from  Section  4  that  the  approximations 
to  both  the  open-loop  semigroup  and  its  adjoint  converge  strongly.  Also, 
recall  that  Hypothesis  8.2  guarantees  a  unique  nonnegative  selfadjoint  solu¬ 
tion  to  the  Riceati  equation  (8.5)  for  each  n.  Replacing  A  and  B  with  A* 

n  n  n 

* 

and  C  in  Theorems  5.1  and  5.3,  we  obtain 
n 

a 

Theorem  8.3.  i)  If  II  1^  1 1  is  bounded  uniformly  in  n,  then  the  Riccati 

A  A 

algebraic  equation  (7.16)  has  a  nonnegative  selfadjoint  solution  II  and  ^En 
converges  weakly  to  ffc.  ii)  If  there  exist  positive  constants  M  and  p, 
independent  of  n,  such  that 


It)  I|  <  Me"pt  ,  t^O, 

(8.28) 

then  ||  1 1  is  bounded  uniformly  in  n,  K  Pgn  converges  strongly  to  II  and 

a  **-i  A 

exp( [An  -W  Cn]t  )PEn  converges  strongly  to  S(t),  the  semigroup  generated 

A  *A_i  A 

by  A- DC  R  C,  the  convergence  uniform  in  t  1  0.  iii)  If  is  bounded  away 

A 

from  zero  uniformly  in  n,  then  IH3^||  being  bounded  uniformly  in  n  guarantees 
the  existence  of  positive  constants  M  and  p  for  which  (8.28)  holds  for  all  n.  O 


I  lexpUA  -  C*R_1C 
n  n  n  n 


< 


•■i 
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The  proof  of  the  following  theorem  is  practically  identical  to  that  of 


Theorem  5 .4 . 


Theorem  8.4.  If  Q  is  E-coercive  and  dQ  =  o,  then  there  is  no  nonnegative  sel- 
fadjoint  solution  of  the  Riccati  operator  equation  (7.16),  and 

II  EL  ||  — >  ®  as  n  — >  ®.  O 

^  (8.29) 

Our  purpose  for  bothering  to  state  this  obvious  dual  result  is  to  point  out 
the  following  question.  Can  Theorem  8.4  be  modified  to  include  the  case  where 
Q  has  the  form  (7.23)  with  QQ  =  0 ,  Q1  =  0  and  Q2  coercive  on  H? 

Next,  we  have  the  dual  result  to  Theorem  5.6: 


Theorem  8.5,  Suppose  that  A^  has  an  invariant  subspace  Vq  which  is  also  invari¬ 
ant  under  the  damping  map  Ay^A^,  that  Eq  =  Vq  X  VQ  is  an  observable  subspace, 
and  that  the  restrictions  of  AQ  and  d0(*,»)  to  are  both  H-coercive.  Also, 

suppose  that  VQ  has  finite  dimension  n0  and  that,  for  each  n  2  nfl  in  the 
approximation  scheme,  the  first  rig  e^s  span  VQ  and  the  rest  are  orthogonal  to 
VQ  in  both  V  and  H. 

/S  A 

i)  Then  (7.16)  has  nonnegative  solution  II,  and  HUnll  is  bounded  uniformly 

in  n,  so  that  converges  to  fi  weakly. 

%  X  A 

ii)  If  Eq  and  EQ  (the  E-orthogonal  complement  of  Eq)  are  invariant  under  Q, 

A 

and  if  the  E^-part  of  the  system  (A,Q)  is  controllable,  then  the  hypothesis  of 
Theorem  8.3  ii)  holds. 


Proof.  The  proof  is  practically  identical  to  that  of  Theorem  5.6  with  B 
replaced  by  C*.  For  ii)  ,  note  that,  when  we  partition  A  and  $  as  in  (3.16), 
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the  finite  dimensional  system  is  controllable  if  and  only  if  the 

system  is  observable.  O 


Remark  8.6.  Remarks  5.7  and  5.8  pertain  to  Theorem  8.5  as  well  as  to  Theorem 
5.6;  i.e. ,  in  most  applications  the  theorem  requires  either  that  both  kQ  and 
dp  be  coercive  (so  that  Vq  =  (0})  or  that  the  natural  mode  shapes  be  the  basis 
vectors  and  the  damping  not  couple  the  natural  modes.  It  seems  unlikely  that 
a  finite  number  of  observable  rigid-body  modes  could  change  the  nature  of  the 
convergence,  but  they  greatly  complicate  the  proofs.  For  applications  where 
both  rigid-body  displacement  and  rigid-body  velocity  are  measured,  a  result 
analogous  to  Theorem  5.9  can  be  obtained,  but  we  will  not  bother  here  because 
it  adds  no  significant  insight  and  we  cannot  use  it  in  the  example  in  Sections 
6  and  10.  Also,  see  Remark  10.1.  □ 


Theorem  1.2. 


/v 

converges  strongly  to  II, 


then 


(8.30a) 


* *®in  “  I Ih  °*  as  m 

A  a  A 

where  ^  and  g.^  are  the  functional  estimator  gains  in  (7.20)  and  fin 
are  the  approximating  functional  gains  in  (8.25). 


(8.30b) 

.  * 

and  «in 


Proof.  The  result  follows  from  (7.22)  and  (8.26).  Q 
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9.  The  Finite  Dimensional  Compensators  and  Realizable  Closed-loop  Systems 


9.1  Closing  the  Loop 


lu  ^  Vi 

The  ntn  compensator  consists  of  the  ntn  approximation  to  the  optimal  con¬ 
trol  law  in  Section  4,  applied  to  the  output  of  the  nth  estimator/observer  in 
Section  8;  i.e. ,  the  feedback  control 


'  Fn*n 


(9.1) 


where 


Pn  '  "'SX 


(recall  (4.25))  and  zR(t)  is  the  solution  to  (8.3) 


(9.2) 


Equivalently,  this  com¬ 


pensator  can  be  written  as 


un  =  Fnn 


(9.3) 


where 


?n  =  R~1Bn'nn 


(9.4) 


(recall  (4.43))  and  the  2n-vector  *(t)  is  the  solution  to  (8.6)).  On-line  com¬ 
putations  will  be  based  upon  the  latter  representation,  and  the  block  diagram 
in  Figure  9.1  shows  the  realizable  closed-loop  system  that  results  from  the 
nth  compensator.  We  will  refer  to  this  system  as  the  £th  closed-loop  system. 
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r 


P 


L 


Note  that  kmr  q  has  compact  resolvent  if  and  only  if  A  does. 

9.2  Convergence  of  the  Closed-loop  Systems 

Now  we  will  consider  the  sense  in  which  the  nfc^  closed-loop  system 
approximates  the  optimal  closed-loop  system  in  Section  7  (Definition  7.3). 
Recall  from  Sections  4.1  and  8.1  how  the  approximating  open-loop  semigroups 

Tn(*)  and  their  adjoints  converge  strongly  and  how  the  input  operators  Bn,  the 
measurement  operators  Cn  and  their  respective  adjoints  converge  in  norm.  Sec¬ 
tions  5  and  8  have  given  sufficient  conditions  for  the  approximating  control 
and  estimator  gains  to  oonverge  to  the  gains  for  the  optimal  infinite  dimen¬ 
sional  compensator.  In  this  section,  we  will  assune 

Hypothesis  9.1 .  As  n  — »  «. 

"F„pEn  *  Fl'~>  “• 

0  (,<H 

Remark  9.2 .  Of  course,  we  are  interested  primarily  in  the  case  where  the 

A  A 

gains  F  and  F  are  the  optimal  LQG  gains  in  (7.15)  and  (7.17)  and  Fq  and  pn  are 
the  corresponding  approximations  in  Sections  4  and  8  (i.e. ,(9.2)  and  (8.4)). 
However,  for  the  analysis  of  this  section,  we  need  only  Hypothesis  9.1  for 
some  F  6  L(E,R°),  F  e  L(rP,E)  and  approximating  sequences  Fr  and  !Fn.  Any  such 
gain  operators  will  yield  closed-loop  semigroup  generators  A^,^  in  (7.9)  and 
A..n  in  (9.6) .  □ 
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We  denote  the  projection  of  E  X  E  onto  EXE  by  P„-  . 

n  *  EEn 

Theorem  9.3.  For  t  2  0,  s«,»n(t)PEEn  oomrerges  strongly  to  ,  and  the 

convergence  Is  uniform  in  t  for  t  in  bounded  intervals. 

Proof.  This  follows  from  the  strong  convergence  of  the  open-loop  semigroups 
and  the  uniform  norm  convergence  of  the  control  and  estimator  gains.  □ 


We  should  expect  at  least  Theorem  9.3,  but  we  need  more.  We  should 
require,  for  example,  that  if  S(t)  is  uniformly  exponentially  stable,  then 
S  .  <t)  must  be  also  for  n  sufficiently  large.  Although  nunerioal  results  for 
numerous  examples  with  various  kinds  of  damping  and  approximations  suggest 
that  this  is  usually  true,  we  have  been  unable  to  prove  it  in  general.  We  do 
have  the  result  for  the  following  important  case. 


Suppose  that  the  basis  vectors  e^  of  the  approximation  scheme  are  the 

natural  modes  of  undamped  free  vibration  and  that  the  structural  damping  does 

not  couple  the  modes.  Then,  for  each  each  n,  E  and  E^reduoe  the  open-loop 

n  n 

semigroup  T(t)  and  its  generator  A.  For  this  case,  we  can  extend  to  D(A)X 

D(A)  as 


a  _BFnPEn] 
”n  "  lFnC  ^nCompJ 


(9.10) 


where 


AnComp  '  lAn'BnFn~FnCn]PEn  +  A*D(A)fl  E1' 

n  (9.11) 

Note  that  E'f'is  the  span  of  the  modes  not  represented  in  the  nth  compensator. 


The  operator  A^,^  generates  a  semigroup  Se9,n(t)  or  E  X  E,  E  X  and 
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{0}XEn  reduce  and  the  restriction  of  S^.^t)  to  E  X  En  is  S^.^t). 

Hence  S^,  (t)  is  uniformly  exponentially  stable  if  and  only  if  both  S<n>n(t) 
and  the  part  of  the  open- loop  system  on  Ej^  are  uniformly  exponentially  stable. 

Theorem  9.4.  i)  Suppose  that  the  basis  vectors  of  the  approximation  scheme 
are  the  natural  modes  of  undamped  free  vibration  and  that  the  structural  damp- 
ing  does  not  oouple  the  modes.  Then  Seo«n(t)  converges  in  norm  to  S^.^ft), 
uniformly  in  bounded  t- intervals. 

ii)  If,  additionally,  S^.^tt)  is  uniformly  exponentially  stable,  then 
is  uniformly  exponentially  stable  for  n  sufficiently  large. 


Proof.  From  (9.6),  (9.10)  and  (9.11),  we  have 


-A 


oo'  n 


.  [F-Fn1C  4„ 


where 


(9.12) 


=  <BnFnPEn-BF) 


(F  C  P„  -FC) . 
n  n  bn 


Therefore,  I  f A  ,  —A  ,  II  — >  0  as  n  — >®,  and  the  theorem  follows.  □ 

CD  QD  CD  II  * 


(9.13) 


This  paper  emphasizes  using  the  convergence  of  the  approximating  control 

A 

and  estimator  gain  operators  Fn  and  Fn,  and  the  convergence  of  the  functional 
gains  that  can  be  used  to  represent  these  operators,  to  determine  the  finite 
dimensional  compensator  that  will  produce  essentially  optimal  closed-loop  per¬ 
formance.  However,  close  examination  of  the  right  sides  of  (9.12)  and  (9.13) 
reveals  another  important  convergence  question.  While  the  gain  convergence  in 
(9.8)  and  (9.9)  drives  the  off-diagonal  blocks  in  (9.12)  to  zero,  the  norm 
convergence  of  the  approximating  input  and  output  operators  also  is  essential 
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in  killing  An>  Expanding  the  two  terms  in  this  block  yields 


BnFnPEn 


-  BF 


=  B_(FnP-  -F) 
n  n  En 


<Bn-B)F, 


(9.14) 


A  A 

F  C  Pc  -  FC 
n  n  En 


(VF)CnPEn 


F(CnPEn-C)‘ 


(9.15) 


The  second  term  on  the  right  side  of  each  of  these  equations  represents, 
respectively,  control  and  observation  spillover,  which  has  been  studied  exten¬ 
sively  by  Balas.  Together,  the  control  spillover  and  observation  spillover 
oouple  the  modes  modelled  in  the  compensator  with  the  modes  not  modelled  in 
the  compensator.  The  spillover  must  go  to  zero  —  as  it  does  when  B^  and  Cn 
converge  --  for  A  -  T  to  go  to  zero. 

i  #  n 


We  should  ask  then  whether  there  exists  a  correlation  between  the  con- 

A 

vergence  of  Fr  and  Ffi  and  the  elimination  of  spillover.  The  answer  is  ye3  if 
no  modes  lie  in  the  null  space  of  the  state  weighting  operator  Q  in  the  per¬ 
formance  index  and  if  the  assumed  process  noise,  whose  covariance  operator  is 

A 

Q,  excites  all  modes,  but  this  correlation  is  difficult  to  quantify.  As  we 
discussed  in  Section  6.4,  the  two  main  factors  that  determine  the  convergence 
rates  of  the  gains  are  the  Q-to-R  ratio  and  the  damping,  neither  of  which 
affects  the  convergence  of  B^  and  C  .  On  the  other  hand,  when  either  factor 
(small  Q/R  or  large  damping)  causes  the  gains  to  converge  fast,  it  generally 

A 

also  causes  the  magnitude  of  F  and  F  to  be  relatively  small,  thereby  reducing 
the  magnitude  of  the  spillover  terms  in  (9.14)  and  (9.15).  Also,  as  n 
increases,  the  increasing  frequencies  of  the  truncated  modes  usually  reduce 
the  ooupling  effect  of  spillover.  This  is  well  known,  although  it  cannot  be 
seen  from  the  equations  here.  In  examples  that  we  have  worked,  we  have  found 

that  when  n  is  large  enough  to  produce  convergence  of  the  control  and 
estimator  gains,  the  effect  of  any  remaining  spillover  is  negligible.  But 

this  may  not  always  be  true,  and  spillover  should  be  remembered. 
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9.3  Convergence  of  the  Compensator  Transfer  Functions 


The  transfer  function  of  the  nth  compensator  (shown  in  the  bottom  block 
of  Figure  9.1)  is 


fi 

n 


(s) 


"  -Fn<sI  ‘  tAn-BnFn 


+F, 


(C 


0Fn~Cn)I} 


n' 


(9.16) 


which  is  an  m  X  p  matrix  function  of  the  oomplex  variable  s  for  each  n,  as  is 
the  similar  transfer  function  ff(s)  in  (7.14)  for  the  infinite  dimensional  com¬ 
pensator.  We  continue  to  assume  Hypothesis  9.1. 


We  will  denote  the  resolvent  set  of 
P([A-BF+F(C0f-C)]) . 


£A-BF+$(CnF-C)] 


by 


Theorem  9.3.  There  exists  a  real  number  Sj  such  that,  if  Re(s)  >  a1#  then  s  e 
p(IAn~Bn  Fn+Fn(C0Fn~Cn)] )  for  all  n,  and  ®n(s)  converges  to  Sts),  uniformly  in 
compact  subsets  of  such  s. 

-Ero.of •  The  operator  [A-BF+F (C^f-C)]  is  obtained  from  a  contraction  semigroup 
generator  by  perturbation  with  bounded  operators,  and  the  approximations  to 
the  perturbing  operators  are  bounded  in  n,  by  strong  convergence.  In  view  of 
this,  close  examination  of  the  basic  approximation  scheme  in  Section  4.1  will 
show  that  there  exists  a  bound  of  the  form  M1exp(a1t),  independent  of  n,  for 
the  semigroups  generated  by  Un-Bn  Fn+Fn < C0Fn~Cn) ) .  Also,  these  semigroups 
converge  strongly  to  the  semigroup  generated  by  [A-BF+£(C0f-C)J ,  according  to 
[G3,  Theorem  6.6],  For  Re(s)  >  then,  the  resolvent  operator  in  $n(s)  con¬ 
verges  strongly  to  that  in  S(s),  uniformly  in  oompact  s-subsets,  by  IK1,  page 
504,  Theorem  2.16,  and  page  427,  Theorem  1.2].  □ 
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This  result  leaves  much  to  be  desired.  For  example,  it  does  not  guaran¬ 
tee  that  any  subset  of  the  imaginary  axis  will  lie  in 

p^^n-pnpn+pn^0pn-  Cn)])  for  sufficiently  large  n,  even  if  all  of  the  ima¬ 
ginary  axis  lies  in  p([A-BF+F(CqF-C)3 ) .  As  with  the  convergence  of  the 
closed-loop  systems,  we  can  get  more  for  certain  important  cases. 


Remark  9.6.  If  the  open-loop  semigroup  T(*)  (whose  generator  is  A)  is  an  ana¬ 
lytic  semigroup,  then  there  exist  real  numbers  a,  0  and  M,  with  0  and  M  posi¬ 
tive,  such  that  p( [A-BF+£(CqF~c)] )  contains  the  sector  {s : larg(s-a) I <  2  +  » 

and  for  each  s  in  this  sector. 


1 1  (sI-[A-BF+F(C0F-C)])_1||  <  M/|3-a| 


(9.17) 


Theorem  9.1.  i)  If  the  basis  vectors  of  the  approximation  scheme  are  the 
natural  modes  of  undamped  free  vibration  and  the  structural  damping  does  not 
couple  the  modes,  then  each  3  in  p( [A-BF+F(C0f-C>] )  is  in 

p( [An-BnFn+Fn(CoFn-Cn)] )  for  n  sufficiently  large  and  Bn(s )  converges  to  Sts) 

A 

as  n  — >  «  ,  uniformly  In  compact  subsets  of  p(lA-BF+F(C0F-C)l ) .  ii)  If, 
additionally,  T(«)  is  an  analytic  semigroup,  then  ®n(s)  converges  to  S(s)  uni¬ 
formly  in  the  sector  described  in  Remark  9.6. 


Proof,  i)  In  this  case,  we  have  also 


*»<•>  - 


FnpEn^sI  ”  ^nComp^  F 


n' 


(9.18) 


where  ^ncomp  is  the  operator  on  D(A)  defined  by  (9.11).  The  result  follows 
from  (9.8)  and  (9.9)  and  the  fact  that  ^nQomp  converges  in  norm  to 

r  A 

IA-BF+F(Cqf-C)] .  ii)  The  result  follows  from  i)  and  a  bound  on 

.  mi  _1 

(sI-Ancomp)  for  large  |s|  that  is  obtained  from  the  Neumann  series  in  view 
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of  (9.17)  and  the  uniform-norm  convergence  of 

Theorem  9.8.  If  A  has  oompact  resolvent,  then  ®n(s)  converges  to  8(s)  for 
each  s  e  p([A-BF+F(CgF-C)] ) ,  uniformly  in  oompact  subsets. 

Proof.  As  a  result  of  Theorem  4.4,  the  resolvent  operator  in  *n(s)  converges 
in  norm  to  the  resolvent  operator  in  S(s)  for  sufficiently  large  real  s. 

After  an  artificial  extension  of  Ar  to  En  then,  the  present  theorem  follows 
from  [Kl,  pages  206-207,  Theorem  2.25]. 
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10.  Closing  the  Loop  in  the  Example 

As  in  Definition  7.3,  the  optimal  closed-loop  system  is  formed  with  the 
optimal  infinite  dimensional  compensator,  which  consists  of  the  optimal  con¬ 
trol  law  for  the  distributed  model  of  the  structure  applied  to  the  output  of 

an  optimal  infinite  dimensional  state  estimator.  This  optimal  control  law  is 
the  limit  of  the  approximating  finite  dimensional  control  laws  in  Section  6. 

In  this  section,  we  first  approximate  the  infinite  dimensional  estimator,  as 

in  Section  8,  and  then  apply  the  approximating  control  laws  in  Section  6  to 
the  approximating  finite  dimensional  estimators  to  produce  a  sequence  of  fin¬ 
ite  dimensional  compensators  that  approximate  the  optimal  compensator. 


10.1  The  Estimator  Problem 


We  assume  that  the  only  measurement  is  the  rigid-body  angle  0  and  that 
this  measurement  has  zero-mean  Gaussian  white  noise  with  variance  $  =  10~4. 

We  model  the  process  noise  as  a  zero-mean  Gaussian  white  disturbance  that  has 
a  component  distributed  uniformly  over  the  beam,  as  well  as  two  concentrated 
components  that  exert  a  force  on  the  tip  mass  and  a  moment  on  the  hub.  For 

A  A 

this  disturbance,  the  covariance  operator  Q  has  the  form  (7.23)  with  QQ=  o, 

\  =  0  and  Q2  =  I. 

We  construct  the  approximating  estimators  as  in  Section  8.1.  The  gain 
for  the  nth  estimator  is  given  by  (8.13)  with  the  solution  to  the  Riccati 
matrix  equation  (8.14).  For  the  rigid-body  measurement,  the  matrix  Cn  is 


Cn  =  [1  0  0  0...]. 


(10.1) 
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According  to  (8.15),  the  matrix  Qn  is 


(10.2) 


since  Wn  is  the  matrix  in  (4.33).  (As  always,  *Cn  is  the  inverse  of  the  mass 

matrix.)  Recall  from  Section  6.3  that  n  =  2n  +  l  where  n  is  the  number  of 

©  © 

elements. 


Our  only  use  for  the  functional  estimator  gains  is  to  measure  the  conver¬ 
gence  of  the  finite  dimensional  estimators  to  the  optimal  infinite  dimensional 
estimator.  To  see  the  convergence  of  the  approximating  estimator  gains,  we 
compute  the  approximating  functional  estimator  gains  as  in  Section  8.3.  Like 
the  functional  control  gains,  the  functional  estimator  gains  have  the  form 


^  * 

«  "  (V*g,Pg)  ' 

and  the  corresponding  approximations  have  the  form 


(10.3a) 


(10.3b) 


A 

f 


'  (afn,0gn’Pgn)  * 


(10.4a) 


A 


^“gn’^gn’Pgn*  * 


(10.4b) 


Remark  10.1  We  cannot  guarantee  as  much  about  convergence  for  the  approximat¬ 
ing  estimators  as  we  could  for  the  approximating  control  problems  in  Section 
6.  Since  the  damping  in  this  example  does  not  couple  the  natural  modes  and 
the  rigid-body  mode  is  observable,  we  would  have  Part  i)  of  Theorem  8.5  if  we 
were  using  the  natural  mode  shapes  as  basis  vectors.  Therefore,  we  know  at 
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least  that  a  solution  to  the  Infinite  dimensional  Riccati  equation  (7.16) 
exists  and  that  the  infinite  dimensional  estimator  that  we  want  to  approximate 
exists.  The  mmerical  results  indicate  that  the  solutions  to  the  finite 
dimensional  Riccati  equations  are  bounded  in  n  and  that  the  functional  esti¬ 
mator  gains  converge  in  norm.  The  rigid-body  mode  prevents  our  guaranteeing  a 
priori  all  the  convergence  that  we  want.  If  a  torsional  spring  and  damper 
were  attached  to  the  hub  in  the  current  example,  we  would  have  ooercive  stiff¬ 
ness  and  damping  and  Theorem  8.5  il)  would  guarantee  that  the  solutions  to  the 
finite  dimensional  Riccati  equations  converge  strongly  and  that  the  functional 
estimator  gains  converge  in  norm  for  the  basis  vectors  used  here.  Also,  see 
Remark  6.2  and  Remark  8.6.  □ 

—  A  9  • 

For  damping  coefficient  cQ  =  10  ,  Figures  10.1  and  10.2  show  dfn  and  and 

Tables  10.1  and  10.2  list  the  the  scalars  a-  ,a  and  .  Since  d*>  (0)  = 

rn  gn  'gn  *in 

d^n<0)  =  0,  the  convergence  of  df„  implies  the  convergence  of  pfn  = 

as  in  the  control  problem,  Pfn  is  not  an  independent  piece  of  information 
about  the  estimator  gains  while,  as  far  as  our  results  go,  8  is.  We  main- 
tain  analogy  with  the  control  problem  and  list  only  Pgn  in  the  table. 


« 


« 
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number  of  elements  n  =  2, 3, 4, 5, 8  number  of  elsnent 


number  of  elements  n,  =  4 .  6,  8,  10  number  of  elements 


ne 

°fn 

agn 

^gn 

2 

5.0358 

12.680 

-1334.9 

3 

5.2514 

13.789 

-1455.4 

4 

5.3195 

14.149 

-1495.7 

5 

5.3478 

14.300 

-1512.2 

8 

5.3611 

14.371 

-1520.1 

Table  10.1.  Scalar  Components  of  Functional  Estimator  Gains 

—4  —4 

Damping  coefficient  c^  =  io  ;  estimator  R  =  10 
number  of  elements  ng  =  2,  3,  4,  5,  8 


°fn 


gn 


4  5.3195  14.149 

6  5.3567  14.347 

8  5.3611  14.371 

10  5.3623  14.377 


-1495.7 

-1517.5 

-1520.1 

-1520.8 


Table  10.2.  Scalar  Components  of  Functional  Estimator  Gains 


Damping  coefficient  cQ 


estimator  R  =  10 


-4 


number  of  elements  ng  =  4,  6,  8,  10 
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Figures  10.3a  and  10.3b  and  Table  10.3  give  the  numerical  results  for  the  fin¬ 
ite  dimensional  estimators  when  the  structural  damping  is  zero.  While  Theorem 
8.4  says  that  the  solutions  to  the  finite  dimensional  Riceati  equations  for 

A 

these  estimators  will  not  converge  when  the  damping  is  zero  and  Q  is  coercive 
on  E,  we  have  no  result  to  predict  the  convergence  for  zero  damping  when  Q  is 

not  coercive  (even  though  is  ooercive  on  H) .  From  the  numerical  results 

A 

though,  f  does  not  appear  to  converge. 


ne 

afn 

a__ 

gn 

^gn 

2 

5.0730 

12.868 

-1354.4 

3 

5.3390 

14.253 

-1506.0 

4 

5.4417 

14.806 

-1568.0 

5 

5.4  894 

15.067 

-1596.3 

8 

5  .5398 

15.345 

-1627.2 

r 

w 

Table  10.3 .  Scalar  Components  of  Functional  Estimator  Gains 
Zero  damping;  estimator  R  =  10-4 

r  t 

number  ng  =  2,  3,  4.  5,  8 


1 
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10.2  Approximation  of  the  Optimal  Compensator 

—4 

Finally,  for  the  damping  Cq  =  10  .  R  =  *05  in  the  control  problem  and  R 

-4 

=  10  in  the  estimator  problem,  we  construct  the  finite  dimensional  compensa¬ 
tor  in  Figure  9.1;  i.e,  for  each  n  =  2ng  +  1,  we  apply  the  n^h  control  law 
represented  by  the  functional  gains  in  Figure  6.4  and  Table  6.4  to  the  output 
of  the  n^h  estimator  represented  by  the  functional  gains  in  Figure  10.2  and 
Table  10.2.  As  the  number  of  elements  increases,  the  transfer  function  in 
(9.16)  of  the  finite  dimensional  compensator  converges  to  the  transfer  func¬ 
tion  in  (7.14)  of  the  optimal  infinite  dimensional  compensator,  as  described 
in  Section  9.3.  Theorem  9.5  and  Remark  9.6  apply.  Figure  10.4  shows  the  fre¬ 
quency  response  (bode  plots)  of  the  finite  dimensional  compensators  for  4,  6, 

8  and  10  elements.  The  phase  plot  is  for  10  elements  only.  These  plots  indi¬ 
cate  that  the  finite  dimensional  compensator  for  eight  or  more  elements  is 
virtually  identical  to  the  optimal  infinite  dimensional  compensator,  as 
predicted  by  the  functional  gain  convergence  in  Figures  6.4  and  10.2. 
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number  of  elements 


10.3  Comments  on  the  Structure  and  Dimension  of  the 
Implementable  Compensators 

Though  this  paper  does  not  address  the  problem  of  obtaining  the  lowest- 
order  compensator  that  closely  approximates  the  Infinite  dimensional  compensa¬ 
tor,  we  should  note  that  the  compensators  based  on  eight  and  ten  elements  here 
are  unnecessarily  large  because  the  finite  element  scheme  that  we  chose  is  not 
nearly  the  most  efficient  in  terms  of  the  dimension  required  for  convergence. 
(The  dimension  of  the  first-order  differential  equation  in  the  compensator  is 
2(2ne+l).)  We  used  cubic  Hermite  splines  here  to  demonstrate  that  the  finite 
element  scheme  most  often  used  to  approximate  beams  in  other  engineering 
applications  can  be  used  in  approximating  the  optimal  compensator.  In  [GS]  , 
we  compare  the  present  scheme  with  one  using  cubic  B-splines  and  one  using  the 

natural  mode  shapes  as  basis  vectors.  The  natural  mode  shapes  yield  the 

fastest  converging  compensators,  but  the  B-splines  are  almost  as  good.  The 

only  advantages  of  the  Hermite  splines  result  from  the  fact  that  the  coding  to 

build  the  basic  matrices  (mass,  stiffness,  etc.)  is  simpler  than  for  B-splines 
and  the  fact  that,  before  the  Riccati  equations  based  on,  say,  ten  natural 
modes  are  solved,  a  much  larger  finite  element  approximation  of  the  structure 
must  be  used  to  get  the  ten  modes  accurately. 

To  understand  the  redundancy  in  the  large  finite  dimensional  compensators 
here,  it  helps  to  consider  the  structure  of  the  optimal  compensator.  It  is 
based  on  an  infinite  dimensional  state  estimator  that  has  a  representation  of 
each  of  the  structure's  modes.  In  the  present  example,  the  optimal  compensa¬ 
tor  estimates  and  controls  the  the  first  six  modes  significantly,  the  next 
three  modes  slightly,  and  virtually  ignores  the  rest.  This  observation  is 
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based  on  the  projections  of  the  functional  gains  onto  the  natural  modes  and  on 
comparison  of  the  open-loop  and  closed-loop  eigenvalues.  (See  [G5]  for  more 
detail,  including  the  spectrum  of  closed-loop  system  —  which  is  stable  — 
obtained  with  the  ten-element  compensator  here.)  The  infinite  dimensional 
compensator  then  has  an  infinite  number  of  modes  that  contribute  nothing  to 
the  input-output  map  of  the  compensator.  These  inactive  modes  are  just  copies 
of  all  the  open-loop  modes  past  the  first  nine.  They  can  be  truncated  from 
the  compensator  without  affecting  the  closed-loop  system  response  signifi¬ 
cantly.  The  nunber  of  active  modes  in  the  compensator  —  i.e.  ,  the  modes  that 
contribute  to  the  input-output  map  —  depends  on  the  structural  damping  and 
the  Q’s  and  R's  in  the  LQG  problem  statement.  (See  the  discussion  in  Section 
6.4  about  the  effect  of  damping  and  control  weighting  on  performance.) 

The  compensator  computed  here  based  on  ten  elements  has  21  modes 
(although  we  did  not  do  the  computations  in  modal  coordinates).  Nine  of  these 
compensator  modes  are  virtually  Identical  to  the  nine  active  modes  in  the 
infinite  dimensional  compensator,  and  the  twelve  inactive  modes  are  approxima¬ 
tions  to  the  tenth  through  twenty-first  open-loop  modes  of  the  structure.  The 
inactive  modes  result  from  the  large  number  of  elements  needed  to  approximate 
the  active  compensator  modes  accurately.  Now  that  we  essentially  have  the 
optimal  compensator  in  the  ten-  element  compensator,  we  could  truncate  the 
twelve  inactive  modes  and  implement  a  compensator  with  nine  modes.  And  we 
probably  could  reduce  the  compensator  even  further  using  an  order  reduction 
method  like  balanced  realizations. 
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11.  Conclusions 


For  the  deterministic  linear-quadratic  optimal  regulator  problem  for  a 
flexible  structure  with  bounded  input  operator  (the  j.n  (2.1)),  the  approxi¬ 
mation  theory  in  Sections  4  and  5  is  reasonably  complete.  The  most  important 
extensions  should  be  to  the  corresponding  (very  difficult)  problem  with 
unbounded  input  operator,  for  which  there  exists  little  approximation  theory. 
Because  of  the  different  kinds  of  boundary  input  operators,  stiffness  opera¬ 
tors  and  structural  damping,  all  of  which  must  be  considered  in  detail  when  BQ 
is  unbounded,  it  seems  unlikely  that  the  approximation  theory  for  the 
unbounded- input  case  can  be  made  as  complete  as  the  theory  here. 

The  convergence  results  in  Section  8  for  the  estimation  problem  are  less 
complete  than  those  for  the  control  problem  because  rigid-body  modes  present 
more  technical  difficulties  for  the  proofs  in  the  estimator  case.  However, 
our  analysis  and  numerical  experience  suggest  that  the  difficulties  only  make 
the  proofs  harder  and  that  the  convergence  in  the  estimation  problem  is  ident¬ 
ical  to  the  convergence  in  the  control  problem,  and  that  controllable  and 
observable  rigid-body  modes  make  no  qualitative  difference  in  either  problem. 

Where  we  would  most  like  substantial  improvement  over  the  results  of  this 
paper  is  in  Section  9.2,  which  considers  how  the  approximating  closed-loop 
systems  obtained  by  controlling  the  distributed  model  of  the  structure  with 
the  finite  dimensional  compensators  converge  to  the  optimal  closed-loop  sys¬ 
tem,  obtained  with  the  infinite  dimensional  compensator.  Theorem  9.4  gives  us 
what  we  want  for  problems  where  the  damping  does  not  couple  the  natural  modes 
of  free  vibration  and  the  natural  mode  shapes  are  the  basis  vectors  for  the 


approximation  scheme.  In  particular,  this  theorem  says  that,  if  the  optimal 
closed-loop  system  is  uniformly  exponentially  stable,  then  so  are  the  approxi¬ 
mating  closed-loop  systems  for  sufficiently  large  order  of  approximation.  We 
have  verified  nuaerically  the  stability  of  the  approximating  close-loop  sys¬ 
tems  for  the  example  in  Sections  6  and  10,  where  the  basis  vectors  are  not  the 
modes.  This  example  and  others  have  made  us  suspect  that  Theorem  9.4  is  true 
when  the  basis  vectors  satisfy  Hypothesis  4.1  only  and  when  the  damping  cou¬ 
ples  the  modes. 

Another  possible  approach  to  analyzing  the  convergence  of  the  approximat¬ 
ing  closed-loop  systems  to  the  optimal  closed-loop  system  is  to  use  the 
input-output  description  in  frequency  domain.  Results  like  those  in  Section 
9.3  are  useful  for  this,  although  for  the  closed-loop  stability  we  want,  we 
probably  need  the  transfer  functions  of  the  finite  dimensional  compensators  to 
converge  more  uniformly  on  the  compensator  resolvent  set  than  we  have  proved 
here.  In  our  example.  Figure  10.4  indicates  that  these  transfer  functions 
converge  uniformly  on  the  imaginary  axis,  but  we  have  no  theorem  that  guaran¬ 
tees  this. 
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APPENDIX 


Errata  for  [Gl] 

In  the  first  paragraph  of  the  proof  of  Theorem  2.1  on  page 
689  of  [Gl],  the  first  sentence  should  be: 

If  a  dissipiive  operator  is  invertible,  its  inverse  is 
dissipative . 

At  the  beginning  of  the  fifth  line  of  the  same  paragraph,  the 
expression  (ax+y)  should  be  deleted  the  first  time  it  occurs. 

The  next-to-last  sentence  of  the  paragraph  should  be : 

Hence,  if  a  densely  defined  maximal  dissipative  operator  has 
dense  range,  its  inverse  is  maximal  dissipative. 

The  theorem  is  correct  as  stated. 

In  the  current  paper,  we  use  Theorem  2.1  of  [Gl]  to  conclude 
that  the  operator  A  defined  in  Section  2  is  maximal  dissipative  (see 

O 

(2 . 10) -(2 .12))  and  that  the  operator  A  in  (2.16)  has  a  unique  maximal 


dissipative  extension. 
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